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Chapter  1 


Introduction:  the  project  objectives,  tasks  and  the  report  structure 

Key  words:  optical  vortices,  phase  singularities,  orbital  angular  momentum,  computer- 
synthesized  holographic  optical  elements,  low-mode  optical  fibers,  vortex  modes  of  optical 
fibers,  optical-vortex  fiber  sensors. 

A  wide  class  of  effects  associated  with  phase  singularities  in  optical  wave  fields,  as  well  as 
topology  of  wave  fronts,  became  recently  the  subject  of  singular  optics ,  new  chapter  of  modem 
photonics  [1].  Although  phase  is  an  auxiliary  function  in  electromagnetic  field  description,  it  is 
very  useful  as  it  gives  a  perception  of  a  wave  propagation  and  transformation  along  its  path.  An 
important  relevant  conception  is  a  wave  front,  or  a  surface  of  equal  phase,  usually  associated  with 
a  crest  of  a  wave,  where  the  field  strength  attains  its  highest  value.  The  normal  to  the  wave  front, 
or  phase  gradient  direction,  shows  the  Poynting  vector  direction  of  energy  flow,  at  least  in  the 
geometrical  optics  approximation  [2].  The  wave  fronts  follow  each  other  with  spatial  separation 
in  one  wavelength,  and  between  two  neighbor  crests  there  are  two  surfaces  where  field  strength 
attains  zero  value,  and  one  surface  where  field  reaches  minimum  (negative)  value  (trough).  This 
perfect  regular  motion,  being  tme  for  a  plane  wave,  can  be  violated  for  real  waves.  In  brief,  phase 
of  a  wave  can  experience  a  “7i-jump”,  which  corresponds  to  a  step  on  a  half  of  a  wavelength  in  a 
wave  train,  producing  a  phase  defect  of  a  wave  front  along  a  continuous  line  in  space.  For 
instance,  some  physical  reasons  can  lead  to  local  variation  of  a  phase  velocity  across  a  wave  front. 
Resulting  bending  of  the  wave  front  can  lead  to  its  tear,  and  the  phase  becomes  indeterminate,  or 
singular  along  the  tear.  The  necessary  and  sufficient  condition  for  phase  singularity  to  appear  is 
the  vanishing  of  the  field  amplitude.  According  to  the  classification  introduced  by  Nye  and  Berry 
[3],  a  monochromatic  light  wave  can  possess  two  main  types  of  phase  singularities:  screw  wave- 
front  dislocation  and  edge  dislocation,  while  mixed  edge-screw  dislocation  is  most  common 
situation.  Nowadays  the  term  “ optical  vortex ”  (OV)  [4]  became  widely  used,  reflecting  the 
general  feature  of  phase  singularities:  phase  circulation  around  the  dislocation  line.  Therefore,  a 
pure  screw  dislocation  is  an  example  of  a  “longitudinal”  OV,  and  edge  dislocation  is  the 
“transversal”  OV,  with  respect  to  the  wave  propagation  direction. 

There  is  now  theoretically  predicted  [5]  and  experimentally  verified  [6]  conclusion  on  the 
existence  of  an  orbital  angular  momentum  for  the  longitudinal  OV  beam.  It  could  be  treated  as 
well  as  classical  effect,  caused  by  the  tangential  (azimuth)  component  of  the  Poynting  vector,  as 
inherent  quantum  property  of  a  photon  [7].  The  deep  insight  to  the  problem  seems  to  be 
important,  as  both  explanations  have  some  obscurities. 

Most  important  for  fundamentals  and  applications  are  longitudinal  OVs,  which  possess  zero- 
amplitude  axis  and  helicoidal  wave  front.  A  vortex  core  can  be  described  with  a  factor 
r|m|exp(/rn^),  where  r  is  the  radius,  cp  is  the  azimuth  angle,  and  m  is  the  signed  integer  called  OV 
topological  charge.  Therefore,  a  laser  beam  with  an  axial  OV  looks  like  dark-core  “hollow”  self¬ 
similar  beam,  which  preserves  its  structure  up  to  infinity,  i.e.  to  the  far  field.  The  unique 
properties  of  OVs  open  a  variety  of  possibilities  for  new  applications. 

The  investigations  of  fundamental  properties  and  possible  applications  of  OVs  develop 
rapidly  during  last  years  [8,9].  This  is  clearly  seen  from  the  Proceedings  of  SPIE  vol.  3487  [10] 
and  vol.  4403  [11]  of  the  1st  (1997)  and  the  2nd  (2000)  International  Conferences  on  Singular 
Optics  (Optical  Vortices),  organized  by  Institute  of  Physics,  NAS  of  Ukraine,  sponsored  by  many 
organizations,  including  EOARD.  The  essential  account  in  organization  of  these  conferences  was 
done  by  Prof.  A.  V.  Volyar  with  co-workers  from  Tavrida  National  University  (Simferopol, 
Crimea,  Ukraine),  the  leading  group  in  the  area  of  investigations  and  applications  of  optical 
vortices  in  fiber  optics. 

At  the  same  time  many  principal  questions  of  optical  vortices  understanding  and 
assessment  are  still  not  established.  To  make  an  essential  contribution  to  this  problem  was  the 
goal  of  the  Partner  Project  P-05 1 . 


The  project  objectives  were  twofold: 


-  more  deep  understanding  of  OVs  unique  properties,  including  their  topological  and 
quantum  features,  which  are  new  for  classical  optics; 

-  the  assessment  of  possible  new  applications  of  OVs  mainly  in  optical- vortex  fiber  sensors 
with  high  sensitivity  and  multi-channel  optical  tweezers  for  manipulation  with  a  system  of 
micro-particles. 

The  investigations  were  carried  out  in  two  institutions,  which  gained  a  reach  experience 
in  OVs  and  are  leading  groups  in  singular  optics: 


1.  INSTITUTE  OF  PHYSICS,  NATIONAL  ACADEMY  OF  SCIENCES  OF 
UKRAINE 

Prospect  Nauki  46,  Kyiv-39,  03650,  Ukraine, 

Tel:  +(38044)  265  12  20;  Fax  +(38044)  265  15  89:  e-mail:  |iegrivko@marion. iop.kiev.ua 

2.  TAVRIDA  NATIONAL  UNIVERSITY,  SIMFEROPOL,  CRIMEA,  UKRAINE 

Yaltinskaya  4,  Simferopol,  95007,  Crimea,  Ukraine 

Tel:  +(380652)  23  22  80,  Fax:  +(380652)  23  23  10;  e-mail:  l~ector@ccssu.crimea.ua 


Project  Manager 

Soskin  Marat  S.,  correspondent-member  of  National  Academy  of  Sciences  of  Ukraine,  Professor, 
Head  of  Department  of  Optical  Quantum  Electronics,  Institute  of  Physics,  National  Academy  of 
Sciences  of  Ukraine, 

Tel.:  +(38044)  265  55  63;  Fax  +(38044)  265  15  89:  e-mail:  msoskin@iop.kiev.ua  and 


marats@vortex.kiev.ua 


The  following  activities  were  planned  to  perform  the  project  objectives: 

Task  1.  Theoretical  consideration  of  the  linear  and  nonlinear  optical  effects,  which 
embrace  OV  beams.  This  task  was  planned  to  contribute  to  the  optical  vortices  understanding 
(report  and  publications)  and  has  been  performed  by  IP  NASU. 

Task  2.  Experimental  investigations  of  OV  beams  diffraction  and  interference.  This  task 
was  planned  to  contribute  to  the  assessment  of  possible  applications  of  optical  vortices  (report  and 
publications)  and  has  been  performed  by  IP  NASU. 

Task  3.  Theoretical  and  experimental  research  of  OV  use  in  multi-channel  tweezers.  This 
task  was  planned  to  contribute  to  the  assessment  of  optical  vortices  possible  applications  for 
manipulating  with  system  of  micro-particles  (modeling  and  report)  and  has  been  performed  by  IP 
NASU. 


Task  4.  Theoretical  and  experimental  studies  of  optical-vortex  fiber  sensors.  This  task 
was  planned  to  contribute  to  the  proposed  optical-vortex  sensor  assessment  (report,  test  protocol, 
publications)  and  has  been  performed  by  TNU. 

Task  5.  Consideration  and  computer  simulation  of  applicability  of  short-pulse  lasers  in 
OV  fiber  sensors.  This  task  was  planned  to  contribute  to  the  assessment  of  the  proposed  sensor 
application  area  (report)  and  has  been  performed  by  TNU. 

The  final  report  structure  is  built  in  accordance  with  the  described  scope  of  activities  and 
consists  from  three  chapters  with  text,  figures,  summary  and  general  conclusion. 


It  is  seen  that  the  project  objectives  cover  a  wide  scope  of  manifold  specific  effects  in  free 
space  propagation,  optical  fibers,  etc.  To  facilitate  the  reading  and  evaluation  of  this  huge 
material  each  relatively  independent  section  possesses  its  own  numeration  of  formulae,  figures 
and  references. 
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Chapter  2 

Interference  and  diffraction  phenomena,  nonlinear  effects 
in  light  beams  with  optical  vortices 

2.1.  Introduction 

As  it  was  pointed  out  above,  light  beams  with  optical  vortices  and/or  edge  dislocations  of 
wave  front  possess  quite  new  properties,  which  are  absent  in  light  beams  with  “usual”  smooth 
wave  front.  This  is  connected  mainly  with  the  next  new  features  of  light  beams  with  dislocations 
of  wave  front: 

-  helical  wave-front  shape  for  optical  vortices  or  wave-front  phase  steps  for  edge  dislocations, 

non-zero  orbital  angular  momentum  of  OV  light  beams,  which  conserves  at  free  propagation. 

Our  analysis  of  state-of-art  of  OVs  understanding  at  the  beginning  of  the  project  has 
shown  that  some  basic  optical  phenomena  like  interference  and  diffraction  of  OV  light  beams 
were  not  considered  yet  carefully.  Therefore,  we  tried  to  fill  this  gap  during  fulfillment  of  the 
project  for  better  understanding  of  OV  physics,  according  to  the  tasks  1  and  2. 

The  material  is  arranged  in  the  following  way,  natural  from  our  point  of  view: 

(i)  basic  properties  of  phase  singularities  in  combined  light  beams  with  two  mutually 
coherent  Gaussian  beams  and  unusual  transformation  of  OV  beam  orbital  angular  momentum  and 
total  topological  charge  by  cylindrical  lens  in  the  focal  plane  and  its  vicinity  (2.2  -  2.7), 

(ii)  second  harmonic  generation  with  OV  beams  with  fractional  topological  charges  (2.8). 


2.2.  Transversal  Optical  Vortices  in  a  Light  Wave:  Structure  and 

Topological  Reactions 


The  purpose  of  the  present  paragraph  is  to  report  a  study  of  a  structure  of  wave-front  edge 
dislocations  in  an  interference  field  of  two-dimensional  (2D)  Gaussian  beams  and  to  discuss  the 
topological  reactions,  which  occur  in  the  process  of  saddle-vortex  collision.  The  central  question 
is  to  find  the  condition,  which  determines  the  direction  of  an  “optical  vortex”  circulation  around 
the  dislocation  line,  which  is  a  line  of  zero  amplitude  of  the  field.  In  contrast  to  the  configuration 
used  in  [1],  we  now  simplify  the  problem  to  2D  (cylindrical)  waves  and  therefore  can  operate 
within  scalar  wave  field  approximation. 


Fig.  1.  .  (a)  A  schematic  view  of  a  2D  Gaussian  beam  with  the  waist  located  at  z  =  0  position,  (b) 
Destructive  interference  of  two  2D  Gaussian  beams  with  different  waist  parameters.  Resulting 
amplitude  distribution  in  the  waist  is  shown  by  dashed  line.  Negative  amplitude  denotes  the  n 
phase  shift,  (c)  A  wave-front  sequence  near  the  waist  (<Fi  =  1):  there  is  a  lack  of  one  wave  front 
between  the  dislocation  lines.  The  wave-front  ledge  before  the  waist  is  directed  along  the 
direction  of  the  wave  propagation  and  opposite  after  the  waist. 

The  geometry  of  waves  is  explained  in  Fig.  1.  A  structure  of  a  single  2D  Gaussian  beam 
is  shown  in  Fig.  la,  with  the  beam  waist  at  z  =  0  position.  The  beam  is  polarized  along  7-axis  and 
propagates  along  Z-axis,  with  diffraction  spread-out  in  ±X-axis  directions.  For  this  2D  beam  of 
Gaussian  shape  the  paraxial  solution  satisfying  paraxial  wave  equation  looks  as  follows: 


E(x,z)  =  Eg 


(1) 


kx2 

2  R(z) 


1  z  .  , 

—  arctan - f  +  <P 

2  z* 


\ 
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where  EG  is  the  amplitude  parameter,  w0  is  the  waist  parameter,  k  is  the  wavenumber,  zR  is  the 
Rayleigh  range,  zR  —  kw\  /  2  ,  R{z)  is  the  radius  of  the  wave  front  curvature  in  the  XZ  plane,  w  = 
w0(l+z2/zR2)1/2,  d>  is  a  constant  phase  term. 

The  solution  (1)  is  very  similar  to  the  description  of  3D  Gaussian  beam  (see  2.5),  except 
another  degree  of  amplitude  diminishing  along  the  propagation  path  z,  (w0/w)1/2,  instead  of  (w0/w) 
for  3D  Gaussian  beam,  and  coefficient  Vi  before  the  Gouy  phase  shift  arctan(z/zR). 

To  produce  a  pair  of  edge  dislocations  parallel  to  the  7-axis  in  the  plane  z  =  0,  let  us  make 
a  superposition  of  two  2D  Gaussian  beams  with  different  amplitudes  E\  >  E2 ,  waist  parameters  w2 
>  wi,  and  a  phase  shift  n  between  their  phases  <Fi  and  <E>2,  very  similar  to  how  it  has  been  done  for 
3D  Gaussian  beams  [1]: 


E(x,z  -  0)=  Ex  exp 

(  2  > 
X 

2 

expz'c&j  +  E2  exp 

(  2  \ 
X 

2 

^  1  ) 

l  ^  J 

expz<f>2. 


(2) 


At  the  waist  plane,  both  waves  have  plane  wave  front,  and  their  relative  phases  are  locked  as  <E>i  = 
d>2  +  71,  thus  producing  destructive  interference,  as  shown  in  Fig.  lb. 

The  edge  dislocation  position  (zero  amplitude  of  the  field)  is  determined  by 


*0  =±. 


Me,/e2). 


(3) 


For  the  analysis  below,  we  define  the  ratio  E2!E\  as  the  governing  parameter  77.  A  family 
of  wave  fronts  (crests)  is  shown  in  figurelb  for  77  =  0.83,  w\  =  10,  w2  =  100,  <E>i  =  1.  (Here  and 
below  all  ranges  are  normalized  for  k  =  1).  First,  it  is  seen  from  Fig.  lc  that  there  is  a  lack  of  one 
wave  front  inside  the  interval  bordered  by  edge  dislocations.  Before  the  waist,  the  wave  front  has 
a  ledge  along  the  direction  of  propagation,  after  the  waist  the  ledge  is  directed  opposite.  The  ledge 
is  smoothed  far  from  the  waist  and  the  resulting  lack  of  the  wavefront  results  in  a  small 
wavelength  difference  within  the  interval  between  the  dislocations  and  outside,  or  phase  velocity 
variation.  The  phase  velocity  is  larger  within  the  interval  between  the  dislocations,  the  physical 
reason  is  the  Gouy  effect  [2]:  a  spatially  compressed  light  beam  increases  the  phase  velocity  on 
the  axis.  To  calculate  the  phase  velocity  of  a  2D  Gaussian  beam  in  the  waist,  Eq.  (1)  can  be 
applied: 


v-c\ 


1  +  - 


2  ZB 


2kzr 


(4) 


where  c  is  the  speed  of  light.  There  is  a  distance  xv  where  v  =  c:  xv  =  w0/ V 2  .  When  x  <  xv,  the 
phase  velocity  exceeds  c,  and  vice  versa.  Due  to  this  phase  velocity  variation,  a  wave  front  being 
plane  at  the  waist  gradually  transforms  to  cylindrical  at  the  far  field. 


Fig.  2.  Lines  of  equal  phase  of  the  interference  field  produced  by  two  2D  Gaussian  beams  in  the 
cross-section  XZ  nearby  the  edge  dislocation,  for  77  =  0.83,  <E>i  =  0.  Saddle  point  is  indicated  as  a 
circle,  vortex  as  a  dot.  There  is  a  71-jump  of  the  phase  at  the  zero-amplitude  point  (core  of  the 
vortex). 


Figure  2  shows  a  family  of  equiphase  lines  in  the  cross-section  XZ  in  the  vicinity  of  the 
edge  dislocation,  for  77  =  0.83.  There  is  an  additional  wave  front  outside  the  inner  stripe  between 
the  dislocation  lines,  and  a  point  where  it  enters  the  periodic  sequence:  here,  at  the  phase  saddle, 
lines  of  <F(x,z)  =  0  and  <F(x,z)  =  2tc  join. 

Analysis  shows  that  the  position  of  the  saddle  can  be  as  inside  the  interval  x0,  (xs  <  x0),  as 
outside  it  (xs  >  x0).  The  physical  sense  of  the  saddle  point  is  the  condition  of  the  energy  current 
stagnation:  vortex  flow  is  exactly  compensated  by  the  longitudinal  flow,  and  a  line  of  standing 
wave  is  created.  Therefore  the  vorticity  area  is  ended  in  the  saddle  point,  and  the  whole  torous 
light  current  around  the  zero-amplitude  line  is  localized  in  a  very  narrow  region.  In  comparison 
with  the  wavelength,  the  separation  between  the  vortex  and  the  saddle  is  extremely  small  and 
amounts  of  about  10"2  of  the  wavelength. 
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-2.0  -1.5  -1.0  -0.5  0.0  0.5  1.0  1.5  2.0 


(c) 

Fig.  3.  Collapse  of  edge  dislocations,  (a)  Equiphase  lines  around  the  dislocations,  which  are  close 
to  the  beam  center,  saddle  points  have  joint,  (b)  Vortices  approach  closer,  saddles  are  split  along 
the  Z-axis.  (c)  All  topological  objects  collide  and  annihilate  at  77  =  1. 


With  increasing  77,  x0  diminishes  according  to  Eq.  (3),  and  two  dislocation  lines  come 
closer  as  well  as  saddles.  The  saddles  which  are  located  inside  the  interval  between  the  vortices 
will  meet  at  x  =  0  point  when  77  =  zR2(2kzR\-\)/  zRi(2kzR2-l).  This  is  the  first  stage  of  the  collapse 
(Fig.  3a).  Then  saddles  appear  again  separately,  but  are  located  now  symmetrical  on  the  Z-axis 
(Fig.  3b).  The  maximum  distance  between  the  saddles  will  be  reached  for  77  ~  0.995  and  amounts 
approximately  to  1 .  In  this  moment  the  distance  between  vortices  becomes  equal  to  the  distance 
between  saddles.  With  further  increase  of  77  the  saddles  as  well  as  vortices  come  closer.  Finally, 
vortices  and  saddles  collide  and  annihilate  altogether  at  the  point  x  =  0  (Fig.  3c). 

With  decreasing  of  77,  saddle  point  moves  toward  the  vortex  and  meets  it  when  T|  =  T|0  = 


exp 


2  2 

wx  -  w2 


2(wj2  +  w22) 


(this  condition  corresponds  to  the  coincidence  of  the  dislocation  line  to  the 


line  of  equality  of  phase  velocities  of  the  interfering  waves,  as  will  be  discussed  below).  In  the 
moment  of  this  collision  a  saddle-vortex  point  appears,  which  is  a  “pathological”  topological 
object.  After  the  crossing  event,  the  saddle  appears  outside  the  dislocation  line,  and  thus  the  light 
flow  circulation  around  the  dislocation  reverses  the  direction. 


Fig.  4.  Solid  line:  cross-section  of  the  surface  of  destructive  interference  (71  phase  difference 
between  interfering  beams),  for  x  >  0.  Dashed  line:  cross-section  of  the  surface  of  the  amplitude 
equality.  For  77  >  770,  there  is  only  one  possible  crossing  point  and  therefore  one  zero-amplitude 
line  (edge  dislocation).  For  77  <  770,  there  are  three  points  of  the  crossing  and  therefore  three  edge 
dislocations.  (Saddle  points  are  not  shown  here). 

Let  us  examine  now  how  new  edge  dislocations  can  sprout  from  the  saddle-vortex  line. 
The  cross-section  of  the  surface  of  the  destructive  interference  is  shown  in  Fig.  4  by  solid  line.  To 
find  the  location  of  edge  dislocation  out  of  the  waist  plane,  the  surface  of  the  amplitude  equality 
can  be  found  also,  as  shown  in  Fig.  4  by  dashed  line.  Its  position  depends  on  the  governing 


parameter  77,  whereas  the  line  of  the  destructive  interference  does  not  depend  on  the  amplitude 
ratio  of  the  beams.  This  line  crosses  the  X-axis  exactly  in  the  point,  which  corresponds  to  the 
equality  of  phase  velocities  of  the  interfering  waves.  This  coincidence  is  hardly  occasional,  and 
can  give  a  hint  for  understanding  of  the  vortex  behavior.  The  consequence  can  be  made:  the 
vortex  rotation  is  directed  toward  the  area  of  higher  phase  velocity. 

Summarizing,  we  have  found  how  the  topological  reactions  of  edge  dislocation  collapse 
and  “unfolding”  [3]  occur  in  an  interference  field  of  two  two-dimensional  Gaussian  beams.  The 
results  obtained  are  valid  as  well  for  three-dimensional  beams. 

The  rule  for  determination  of  the  optical  vortex  circulation  direction  around  an  edge 
dislocation  was  established:  the  light  vortex  rotates  toward  the  region  of  the  higher  phase 
velocity.  This  rule  can  be  used  for  the  explanation  of  the  natural  structure  of  phase  circulation 
around  dark  rings  (circular  edge  dislocations)  observed  by  Karman  et  al.  [4].  The  physical  reason 
of  this  feature  is  not  clear  enough  and  this  fundamental  property  of  light  propagation  seems 
worthwhile  for  further  investigations. 
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2.3.  Transformations  of  Laguerre-Gaussian  modes  carrying 
optical  vortex  by  a  cylindrical  lens  and  orbital  angular 
momentum  conservation 

Even  in  free-space  propagation  singular  beams  may  demonstrate  some  unusual 
properties.  Single-charged  optical  phase  singularities  as  morphological  objects  are  robust  with 
respect  to  perturbations.  For  instance,  addition  of  a  small  coherent  background  does  not  destroy 
single-charged  vortex,  but  only  shifts  its  position  to  another  place  where  the  combined  field 
amplitude  has  a  zero  value.  For  vortices  with  multiple  charge  this  operation  will  split  an  initially 
m-charged  vortex  into  \m\  single-charge  vortices  [1].  Intuition  suggests  that  the  total  topological 
charge  would  be  conserved  in  a  beam  propagating  in  free  space  [2,3].  On  the  basis  of  the  present 
experimental  study,  we  demonstrate  what  happens  in  the  process  of  transformation  of  phase 
singularities. 

The  rotation  of  field  around  the  vortex  axis  causes  a  nonzero  value  of  an  orbital  angular 
momentum  of  a  beam  [4,5].  The  origin  of  the  angular  momentum  may  be  easy  explained  from 
simple  consideration.  As  the  wave  front  of  a  singular  beam  has  helicoidal  shape,  the  Poynting 
vector  ^(p^z),  which  is  perpendicular  to  the  wave-front  surface,  has  the  nonzero  tangential 
component  at  each  point.  This  component  equals  5^,9, z)  =  -mS/kp  in  the  paraxial 
approximation.  The  orbital  angular  momentum  density  is  Mz(p,cp,z)  =  (p/c2  )5,(p  (p,(p,z)  in  the  Z- 
axis  direction.  As  the  Poynting  vector  modulus  is  proportional  to  the  square  of  light  wave 
amplitude,  Soc  \  Es( p,(p,z)  | 2,  we  may  obtain  the  expression  for  Mz(p,(p,z)  of  a  singular  beam  in  the 
Z-axis  direction: 
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(1) 


Mz( p,(p,z)  °C-W2  |  £s(p,<p,z) 

In  more  rigorous  approach,  we  need  to  calculate  a  vector  product  to  determine  the  angular 
momentum  as: 


M=£0r  x  (E  x  B),  (2) 

where  80  is  the  permittivity  of  free  space.  The  relation  between  magnetic  and  electrical  field 
components  in  a  monochromatic  electromagnetic  wave  is  given  by 

icoB  =  V  x  E  (3) 

what  follows  from  Maxwell  equations.  To  calculate  the  product  in  eq.  (2)  we  shall  use  the  full 
form  of  the  field  description  and  write 


E  +  E*  B  +  B  * 

M  =£or  x( - x - ) , 


(4) 


where  the  asterisk  denotes  complex  conjugation. 

Taking  into  account  the  non-zero  z  component  of  the  field,  we  come  to  the  time-averaged 
density  of  orbital  angular  momentum  directed  along  Z-axis  [4,5]: 


Mz  =  ^-e0[x(E*  ^-E^-)-y(E* 

2  CO  dy  dy  dx  dx 


(5) 


The  total  orbital  angular  momentum  Lz  of  the  beam  is  an  integral  over  the  beam  cross-section 


Lz  =  j  \M zdxdy  (6) 

— oo  — oo 


The  density  of  orbital  angular  momentum  Mz  of  a  ’’pure"  singular  beam  is  therefore 
proportional  to  the  topological  charge  multiplied  by  beam  intensity  for  any  circularly  symmetric 
amplitude  distribution,  and  the  total  angular  momentum  Lz  to  the  beam  energy,  as  follows  from 
eq.  (5).  The  total  orbital  angular  momentum  Zz,  as  well  as  the  beam  energy  is  conserved  in  a 
medium  without  losses  [4].  The  orbital  angular  momentum  of  a  singular  beam  does  not  depend  on 
the  frequency  of  light,  what  becomes  evident  if  the  beam  energy  is  expressed  as  the  product  of  the 
number  of  photons  on  the  photon  energy,  /uo[5].  For  more  complicated  beam  structure,  the 
conservation  of  orbital  angular  momentum  is  not  so  evident. 

The  process  of  Laguerre-Gaussian  mode  transformation  by  a  cylindrical-lens  mode 
converter  was  considered  theoretically  by  M.V.  Beijersbergen  et  al.  [4].  It  was  shown  that  a  so- 
called  71  converter  reverses  the  sign  of  an  axial  OV,  and  some  considerations  about  the  transfer  of 
an  orbital  angular  momentum  to  the  optical  system  were  presented  [5,6]. 

We  present  experimental  results  for  a  paraxial  OV  beam  with  topological  charge  m  =  1 
passing  through  a  single  cylindrical  lens.  The  beam  was  created  by  diffractic^,  of  a  collimated 
Gaussian  beam  from  He-Ne  laser  output  on  a  computer-  synthesized  grating1:  The  beam  had 
small  divergence  of  order  of  0.3  mrad  (Fig.l).  Thus  created  OV  beam  was  directed  to  a 
cylindrical  lens. 


*  The  structure  of  computer-  synthesized  grating  is  described  in  details  in  the  Section  2.8. 
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Fig.l  Experimental  setup:  1  and  8 
are  the  semitransparent  mirrors; 

2  is  the  computer-synthesized 
hologram;  3  is  the  cylindrical  lens 
with  focal  length  equal  to  120  cm;  4 
and  5  are  the  100%  mirrors,  6  and  7 
lenses  make  the  beam  expander. 


The  cylindrical  lens  used  for  the  experiment  had  the  focal  length  F=  120  cm.  Just  behind 
the  lens  and  before  the  focal  plane  we  observed  a  single  axial  OV,  which  became  more  and  more 
astigmatic  on  the  path  to  the  focal  plane  (Fig.  2a,b,c).  The  interference  patterns  with  the 
additional  collimated  wave  produced  in  the  beam  expander,  shown  in  Fig.  1,  demonstrate  the 
same  charge  of  the  OV  (Fig.  2d,e,f).  The  most  important  transformations  of  the  focused  beam 
occur  in  the  focal  plane  and  behind  it. 

Theoretical  analysis  predicts  the  appearance  of  an  edge  dislocation  in  the  focal  plane  and 
crossing  the  axial  OV  with  the  edge  dislocation.  After  the  crossing,  the  sign  of  the  OV  is  reversed 
[?] 

We  obtained  in  the  experiment  an  astigmatic  OV  in  the  focal  plane  (Fig. 3  a,c).  Behind  the 
focal  plane  we  have  observed  three  OVs,  two  of  them  possess  the  same  sign  as  the  initial  OV 
(Fig. 3b, d),  they  are  shown  by  arrows  in  Fig.  3d  (A  and  C).  The  central  astigmatic  OV  has  the 
opposite  sign  with  respect  the  initial  one  (B  in  Fig.  3d).  The  utmost  OVs  run  out  after  the  focal 
plane  to  the  regions  of  weak  intensity  and  become  non  detectable. 

The  central  OV  exists  in  the  beam  to  the  far  field.  The  foliation  of  the  lateral  OVs  shows 
that  they  belong  to  a  continuous  dislocation  line,  which  is  not  located  in  the  focal  plane,  but  is 
slightly  bent.  Just  on  the  axis,  this  line  makes  the  edge  dislocation.  In  accordance  with  the  theory, 
we  observed  the  inversion  of  the  topological  charge  when  the  axial  OV  intersects  with  the  edge 
dislocation. 

The  distribution  of  the  orbital  angular  momentum  within  the  beam  transforms  after  passing 
the  cylindrical  lens.  The  theoretical  calculations  give  and  important  result:  the  total  orbital  angular 
momentum  of  the  beam  is  conserved,  even  when  the  vortex  reverses  the  sign  of  the  charge  [7,8]. 
This  is  due  to  the  observed  by  us  spatial  rotation  of  the  beam  as  whole.  This  unexpected 
conclusion  shows  new  features  of  OV  beams,  especially  important  for  the  treatment  of  the  photon 
structure  and  some  possible  applications. 

As  a  summary,  our  experimental  analysis  of  light  beams  carrying  optical  vortices  has 
revealed  the  general  properties  of  the  OV  beam  from  the  point  of  view  of  conservation  law  of  the 
orbital  angular  momentum  [7].  We  have  studied  in  details  the  behavior  of  vortices  in  an  OV  beam 
passed  through  cylindrical  lens.  In  brief,  we  obtained  the  following  results: 

1)  At  the  focal  plane  we  obtained  an  astigmatic  OV.  In  our  experiment  we  did  not  obtain  the 
linear  edge  dislocation. 

2)  We  detected  the  existence  of  three  OVs  just  behind  the  focal  plane.  The  utmost  OVs  have  the 
same  signs  as  initial  OV  and  leave  the  beam  after  the  focal  plane  for  the  regions  of  weak 
intensity.  Central  OV  has  opposite  sign  of  the  topological  charge  with  respect  to  the  initial 
one. 

3)  .Instead  of  the  linear  edge  dislocation,  we  detected  a  mixed  screw-edge  dislocation.  To 
explain  this  effect,  we  can  take  into  account  that  the  initial  OV  beam,  created  by  the 
computer-synthesized  hologram,  was  not  a  pure  Laguerre-Gaussian  mode  [1]. 

4)  The  orbital  angular  momentum  distribution  changes  in  the  process  of  beam  propagation:  the 
component  associated  with  the  helicoidal  shape  of  the  wave  front  has  the  opposite  sign  with 
respect  to  the  part  associated  with  the  spatial  rotation  of  a  beam  [9]. 
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Fig. 2.  Intensity  distributions  and  interference  patterns  for  the  OV  beam  focused  by  the  cylindrical 
lens  with  focal  length  F=  120  cm  at  various  distances  D  from  the  focal  plane:  (a)  D=  -55  cm,  (b) 
D  =  -  40  cm,  (c)  D  =  -25  cm.  Corresponding  interference  patterns  (d,e,f)  were  obtained  with  the 
expanded  Gaussian  reference  wave. 


(b)  (d) 


Fig. 3.  Intensity  distributions  and  interference  patterns  for  the  OV  beam  focused  by  the  cylindrical 
lens  with  focal  length  F=  120  cm  at  the  distances  from  the  focal  plane  D  =  0  (a,c)  and  D  =  +3  cm 
(b,d).  Three  OVs  (A,  B  and  C)  are  seen  in  (d):  A  and  C  have  equal  signs  with  respect  to  the  initial 
OV,  but  B  has  the  opposite  sign. 
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2.4.  Mode  Separator  for  a  Beam  with  an  Off-Axis  Optical 

Vortex 

The  goal  of  the  work  was  an  attempt  to  elaborate  some  experimental  technique  for 
separation  of  "doughnut"  and  ’’normal"  photons,  i.e.  vortex  and  vortex-free  components,  from  a 
combined  beam  with  off-axial  OV.  The  idea  of  the  elaborated  method  is  a  use  of  different  Gouy 
phase  shifts  for  the  modes  with  different  (even  and  odd)  modal  indices  [1]. 

The  Laguerre-Gaussian  (LG)  modes  with  a  nonzero  azimuth  coefficient  /  are  beams  with 
/-charged  axial  OV.  The  structure  of  such  beams  is  described  by  the  expression  for  the  LG  mode 
LGpl  (in  cylindrical  coordinates  p,  (p,  z): 


(1) 


where  ELG  is  the  amplitude  parameter,  w0  is  the  parameter  of  the  beam  size  in  the  waist, 

w  =  w0Dz2/z^  g  is  the  mode  transversal  size  on  a  distance  z  from  the  waist, 

R(z )  =  z  b-  z2J  z2  GJ  the  radius  of  curvature  of  a  wave  front  in  the  cross-section  along  the  beam 

axis,  zR  is  the  Rayleigh  range  equal  to  kw02/2 ,  k  is  the  wave  number,  Llp  is  the  associated 

Laguerre  polynomial,  /  is  the  azimuth  index,  p  is  the  radial  index  of  the  mode.  For  /  =£  0,  there  is 
an  axial  OV  with  topological  charge  m  =  /,  because  the  amplitude  on  the  beam  axis  is  equal  to 
zero,  and  the  phase  has  azimuth  dependence  exp(//<p).  The  index  Q  =  2p  +  |/|  +  1  is  termed  as  the 
index  of  a  mode.  This  index  determines  the  accessory  of  a  mode  to  a  family  of  modes  with 
identical  Q.  Physically,  Q  determines  a  value  of  a  small  addition  to  phase  velocity  of  a  mode, 
which  is  connected  with  the  Gouy  phase  shift  arctan  (z/zR)  [2].  For  the  Gaussian  beam  (p,  1  =  0)  Q 
=  1,  and  the  additional  phase  shift  is  equal  to  -n/2  from  the  beam  waist  to  the  far  field  (in 
comparison  with  a  plane  wave).  ’’Doughnut”  mode  LG0l  (p  =  0,  /  =  1)  has  the  modal  index  Q  =  2, 
and  the  relevant  additional  phase  shift  is  equal  to  -n. 

It  happens  frequently  inconvenient  to  realize  the  laser  generation  on  required  LG  mode  in 
the  laboratory  arrangement.  The  simple  way  to  create  wanted  beams  is  a  diffraction  of  a  Gaussian 
beam  by  a  synthesized  hologram  [3].  Thus,  the  first-  order  diffracted  beam  has  in  the  cross- 
section  a  point  with  phase  singularity,  i.e.  OV  with  the  necessary  sign  and  charge.  The  simplest 
way  to  realize  the  combined  beam  in  practice  is  to  shift  the  center  of  the  hologram,  where  the 
fringe  attains  the  branching,  relatively  to  the  incident  Gaussian  beam  axis. 

Let  us  explain  why  this  method  is  able  to  produce  a  combined  beam.  The  synthesized 
hologram  has  a  form  of  thin  amplitude  grating  with  transmission  T  [4]: 

T(p,(p)  =  T0  +  Tl(p)cos(Kpcos(p  +  Mcp) ,  (2) 

where  T0  is  the  average  value  of  the  hologram  transmission,  7j(p)  is  the  contrast  of  fringes,  K  = 
2  Til A,  where  A  is  the  period  of  the  grating,  M  is  the  charge  of  the  phase  singularity.  The  grating 
with  M  =±1  possesses  the  branching  of  the  central  fringe  like  a  ’’fork”.  This  ’’fork”  is  directed 
upwards  for  positive  M,  and  vice  versa  [3].  Let  us  assume  that  the  contrast  of  fringes  Tx(p)  grows 
linearly  with  a  distance  from  the  center  of  the  grating  and  becomes  equal  unity  on  radius  R  (we 
shall  note,  that  the  fringe  branching  point  does  not  coincide  exactly  with  the  center  of  the  grating, 
being  displaced  on  the  small  quantity  MKI2n  [4]).  Then,  the  distribution  of  the  field  in  a  Gaussian 
beam  (w0 «  R)  diffracted  by  the  grating  with  M  =  1,  attains  just  behind  the  grating  the  following 
form  (in  polar  coordinates): 


Ex  <*-^-exp (—^j  +  i(p), 
R  wn 


or  in  Cartesian  coordinates 


x  +  iy  x2+y2 
E i  ~  ^^exp( - — ) 


R 


(3) 


(4) 


what  corresponds  to  the  LG  mode  LG0l.  If  the  grating  center  is  a  little  bit  shifted  relatively  to  the 
center  of  the  readout  beam,  for  example  on  the  value  x0  <  R  along  the  X-axis,  the  corresponding 
distribution  of  the  field  in  the  first  order  of  diffraction  becomes: 


The  expression  (5)  leads  to  superposition  of  the  mode  LG0l  carrying  OV,  and  the  mode 
LG0°  (Gaussian  beam)  as  the  "background”  with  the  contribution  proportional  to  the  x0.  As  a 
result,  the  combined  beam  is  formed.  Such  simplified  approach  allows  us  to  make  the  important 
conclusion  that  the  structure  of  the  field  in  the  first  order  of  diffraction  can  be  presented  as  a 
superposition  of  an  axial  OV  and  a  vortex-free  component. 

The  experimental  setup  is  shown  in  Fig.  1.  The  initial  beam  with  single  OV  was  obtained 
in  the  first  diffraction  order  of  the  appropriate  synthesized  hologram  ( M  =1).  The  center  of  the 
grating  was  shifted  relatively  to  the  center  of  the  incident  Gaussian  beam  (Fig.  2a).  The  resulting 
position  of  the  OV  in  the  diffracted  beam  did  not  coincide  with  the  center  of  the  beam  (Fig.  2b). 
Then  the  beam  was  launched  into  the  “separator”  as  it  is  shown  in  the  Fig.  1. 


Fig.  1.  Experimental  setup.  1  is  the  He-Ne  laser,  2  is  the  computer-  synthesized  hologram,  3  and  4 
are  identical  beam  splitters,  5  and  6  are  100%  reflective  mirrors,  7  and  8  are  two  confocal  lenses. 


(a)  (b) 


Fig.  2.  (a)  The  orientation  of  the  incident  Gaussian  beam  on  the  grating  with  M  =1,  the  center  of 
the  grating  does  not  coincide  with  the  beam  axis,  (b)  Cross-section  view  of  the  diffracted  beam,  the 
point  of  zero  amplitude  is  shifted  from  the  beam  axis. 

In  one  of  arms  of  the  interferometer  formed  by  identical  beam  splitters  3,  4  and  100% 
reflective  mirrors  5,  6  two  confocal  lenses  7  and  8  were  inserted.  Beam  propagation  between  the 


lenses  results  in  the  incursion  of  an  additional  phase  shift  due  to  the  focusing  of  a  beam  and  its 
subsequent  expansion.  The  resulting  phase  shift  equals  to  the  doubled  Gouy  phase  shift,  and  for  a 
mode  with  index  Q  amounts  -Qn.  If,  for  example,  the  optical  path  difference  between  the  beam 
propagating  to  channel  A  without  reflections  and  the  beam  reflected  to  channel  A  from  the  second 
arm  is  equal  to  an  even  number  of  half-waves  for  a  mode  with  even  Q ,  the  entire  beam  will  be 
directed  to  channel  A,  whereas  radiation  in  channel  B  will  be  absent  due  to  interference 
quenching.  For  such  an  adjusting  of  the  interferometer,  the  optical  path  difference  for  a  mode 
with  odd  Q  is  equal  to  the  odd  number  of  half-waves,  and  the  entire  beam  propagates  through 
channel  B. 

The  necessary  adjusting  of  the  interferometer  was  performed  by  reaching  the  maximum 
beam  dumping  in  the  output  B  for  a  Gaussian  beam  (Q=  1)  entering  into  the  separator. 


(c)  (d) 


Fig. 3.  The  intensity  distributions  in  a  cross-section  view  in  the  output  A  (a)  and  B  (b). 
Interferograms  of  the  beams  in  the  output  A  (c)  and  B  (d). 

Because  the  beam  with  single  off-axis  OV  can  be  represented  as  a  superposition  of  modes 
with  even  Q  (vortex  component)  and  odd  Q  (vortex-free  component),  the  separator  divides  vortex 
and  vortex-free  components  into  different  output  channels.  We  note  that  under  reflection  the  sign 
of  OV  changes  to  the  opposite,  but  due  to  even  number  of  reflections  the  beam  in  the  output  A 
carries  OV  of  the  same  sign  as  the  initial  OV  beam.  With  the  proper  adjusting  of  the 
interferometer,  the  expected  patterns  were  observed  in  the  outputs  A  and  B:  axial  OV  in  the 
output  A  and  a  beam  with  a  smooth  wave  front  in  the  output  B  (Fig.  3a, b).  The  beams  were  tested 
using  interference  with  an  additional  plane  wave  (Fig.  3c, d). 

It  is  clearly  seen  that  by  use  of  a  pure  optical  method  we  succeeded  to  acieve  spatial 
separation  of  photons  belonging  to  the  ’’doughnut”  LG0l  mode  and  to  the  Gaussian  beam  ( LG0° 
mode).  More  general,  this  method  is  suitable  for  separation  of  modes  with  odd  indices  Q  and  even 
0. 
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2.5.  Spontaneous  birth  of  optical  vortices  in  a  system  of  non¬ 
coaxial  co-propagating  Gaussian  beams 


Light  beams  carrying  optical  vortices  were  realized  firstly  in  oscillation  of  “doughnut” 
LG0m  modes  in  He-Ne  laser  with  special  axial  cavity  [1].  Nowadays,  light  beams  with  embedded 
OVs  can  be  obtained  easy  by  the  help  of  computer-synthesized  holograms  which  contains  an 
interference  pattern  of  desired  OV  beam  and  a  reference  beam  [2].  The  survey  on  all  main  devices 
for  vortex  synthesis  can  be  found  in  [3].  We  can  consider  all  these  methods  as  “induced”  birth  of 
vortices  in  an  incident  light  beam  with  smooth  wave  front,  typically  Gaussian  beam. 

But  is  it  possible  to  realize  spontaneous  (natural)  birth  of  vortices  due  to  gradual 
topological  changes  of  smooth  wave  front  during  a  light  beam  free  propagation  or  interaction  with 
nonlinear  media?  This  question  of  principal  importance  for  OVs  understanding  was  answered 
both  theoretically  and  experimentally. The  general  physical  approach  to  solve  the  formulated 
problem  is  based  on  the  consideration  of  combined  light  beams  carrying  optical  vortices  [4]. 

The  theory  of  spontaneous  birth  of  vortices  in  the  superposition  of  co-propagating  but 
non-coaxial  Gaussian  beams  has  been  recently  developed  [5].  Here  we  briefly  summarize  the 
main  conclusions  of  this  theory  to  provide  comparison  with  the  experimental  results.  Theoretical 
treatment  is  based  on  the  model  reducing  the  stationary  three-dimensional  problem  of  propagation 
and  interference  of  light  beams  to  its  two-dimensional  analog,  where  only  transversal  components 
of  a  wave  vector  are  taken  into  account.  The  paraxial  propagation  of  Gaussian  beams  is  described 
by  a  parabolic  equation  that  is  similar  to  the  diffusion  equation,  and  diffraction  spreading  of  such 
beams  turns  out  to  be  the  function  of  the  parameter  z/zR,  where  zR  is  the  Rayleigh  range  (i.e.,  the 
extent  of  a  near-field  domain).  Such  approach  permits  to  find  the  dynamics  of  evolution  of 
transversal  structure  of  interfering  beams  due  to  diffraction.  Computations  show  that  the 
asymmetry  of  an  interference  field  of  the  propagating  beams  plays  the  key  role.  Theory  predicts 
(by  computer  modeling)  the  OV  “dipole”  nucleation  due  to  the  fully  destructive  interference  and 
wave  front  disruption  on  some  distance  behind  the  common  waist  plane. 

To  verify  the  theoretical  predictions,  we  have  realized  the  system  of  two  mutually 
coherent  co-propagating  Gaussian  beams  [6].  The  experiment  was  performed  in  the  scheme  of  a 
twin  Mach-Zander  interferometer  (Fig.l).  He-Ne  laser  (A  =  0.6328  pm,  P  =  40  mW)  was  the 
source  of  the  initial  Gaussian  beam.  The  outer  interferometer  (BS1,  Ml ,  M4,  BS4)  provides 
observation  of  an  interference  pattern  originating  under  superposition  of  the  vortex-bearing  field 
and  expanded  (by  BE)  off-axis  plane  reference  wave.  The  rotation  of  BS4  around  the  X-axis  on 
the  angle  a  permited  to  adjust  the  period  of  interference  patterns.  The  inner  interferometer  ( BS2 \ 
M2,  BS3 ,  and  M3)  creates  two  co-propagating  beams  which  overlap  in  the  near  field.  Both  beams 
after  the  beamsplitter  BS2  undergo  spatial  filtering  by  the  microscopic  lenses  and  pinholes.  In  our 
experiment,  the  pinhole  PHI  was  10  pm  in  diameter,  and  the  pinhole  PH2  was  30  pm  in 
diameter.  The  lenses  LI  (focal  length f_  =  35  mm)  and  L2  (focal  length^  =  50  mm)  transform 
diverging  beams  into  converging  ones  with  the  waists  located  at  the  common  plane  PW.  The 
difference  of  the  focal  lengths  of  the  lenses  LI  and  L2  and  their  distances  from  the  corresponding 
pinholes  (d\  =  37  mm,  d2  =  65  mm)  was  used  to  form  two  coaxial  beams  with  different  waist 
radii.  Displacement  of  the  mirror  M3  along  7-axis  permits  to  shift  the  second  beam  axis  at  any 
distance  A y  .  For  that,  the  beams  leave  to  be  co-propagating. 


Fig.l.  Experimental  arrangement  for  the  observation  of  interference  between  two  co-propagating 
but  non-coaxial  Gaussian  beams:  BS  are  beam  splitters,  Ml-4  are  the  mirrors,  PHI, 2  are  the 
pinholes,  BE  is  the  beam  expander,  LI, 2  are  the  lenses,  PW  is  the  common  waist  plane,  OP  is  the 
observation  plane  where  CCD  camera  was  placed. 

The  measured  spatial  intensity  distributions  at  the  waist  plane  for  the  each  of  two  beams 
and  at  the  plane  distanced  from  the  waist  plane  by  z  =  80  mm  are  shown  in  Fig. 2a  and  Fig. 2b, 
respectively. 


y.  ju  m  y,  jum 


(a)  (b) 

Fig.  2.  .  Intensity  distributions  for  two  overlapping  co-propagating  Gaussian  beams  (•  -  beam  1,  ■ 
-  beam  2):  (a)  at  the  common  waist  plane  (r0 2  /  r0i  =  0.6,  Ay  =  0.4  r0 1);  (b)  the  same  at  the  plane 
distanced  from  the  waist  plane  by  z  =  9.7  zR1. 

Experimental  data  were  fitted  by  Gaussian  envelope  with  good  agreement.  The  waist 
radii  measured  on  the  level  e  ~2  of  the  central  intensity  distribution  maxima  were  r0 1  =83  ±  2  pm 
and  r02  =  0.6  r0 1  =  50  ±  2  pm,  respectively.  The  transversal  beam  displacement  amounted  to  Ay  = 
0.4  r0 1.  The  beam  radii  measured  at  the  level  e  ~2  of  the  central  intensity  distribution  maxima  at  the 
distance  z  =  80  mm  were  rx  =139  ±  2  pm  and  r2  =1 15  ±  2  pm  respectively. 

Using  the  measured  r0 1  and  r02  values,  we  calculated  the  corresponding  Rayleigh  ranges: 
zR\  =  35±2  mm,  and  zR2  =  0.36  zR2  =  12.5  ±1  mm.  The  divergence  of  the  beams  was  measured  to 
be  0i  =  2.4  mrad  and  02  =  4  mrad,  respectively.  The  intensity  maxima  at  the  plane  z  =  80  mm 
were  distanced  by  Ay  =  0.4r0i  =  34±5  pm. 

It  would  be  of  interest  to  observe  in  experiment  the  OVs  trajectory  and  their 
disappearance.  For  this  purpose,  we  obtain  the  set  of  interference  patterns  between  the  combined 
beam  consisting  from  two  Gaussian  beams  and  a  tilted  plane  reference  wave  at  various  distances 
z  from  the  waist  plane  (Fig.  3). 


(d)  (e)  (f) 

Fig.  3.  Interference  between  the  combined  beam  and  tilted  reference  wave,  at  various  distances: 
z  =  0.29Ljh  ( a),  z  =  0.5$Ljh  (b),  z  =  l.l6L ^  (c),  z  =  232L ^  (<7),  z  =  4.08Y^  (e), 

z  =  5.8 Lfti  (f).  Ellipses  show  the  forks  visualizing  the  OVs. 


It  is  seen  in  Fig.  3  that  the  vortex  dipole  is  bom  already  at  the  near  field  (fragment  (a)), 
for  z  =  0.29 Lfti  and  exists  over  rather  large  range  of  distances  (up  to  z  =  4.08Z/^  at  fragment 

(e)).  The  pictures  shown  in  Fig.3  are  obtained  with  step-by-step  change  of  the  scale,  in 
correspondence  with  diffraction  spread-out  of  the  Gaussian  beams.  As  the  distance  z  increases, 
the  vortices  become  more  anisotropic.  At  some  distance  (in  our  experiment,  z  =  5.8Z/^)  the 

vortex  dipole  is  converted  into  the  stmcture  (Fig.3f),  which  reminds  circular  edge  dislocation.  In 
the  calculations,  vortices  annihilate  and  do  not  exist  in  the  far  field. 

To  the  conclusion,  we  have  shown  both  theoretically  and  experimentally  the  possibility  of 
spontaneous  (natural)  birth  of  optical  vortices  in  the  system  of  two-co-propagating  in  free  space 
non-coaxial  Gaussian  beams  with  different  parameters. 
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2.6.  Launching  of  an  Optical  Vortex  Beam  into  a  Leaky 

Planar  Waveguide 


Common  experiments  with  plane  reference  wave,  which  are  good  enough  to  demonstrate 
the  helicity  of  an  OV  beam  wavefront  [1,2],  can  not  give  high  resolution  for  an  accurate 
inspection  of  the  details  such  as  Poynting  vector  orientation. 

The  goal  of  this  communication  is  to  report  the  first  to  our  knowledge  experimental  approach  to 
the  direct  test  of  the  OV  beam  structure  with  a  use  of  a  prism- film  coupler  [3].  This  simple  kind 
of  an  interferometric  device  was  shown  to  be  a  very  sensitive  optical  tool  [4].  The  reason  is  the 
multiple-beam  interference  within  a  planar  waveguide,  producing  high  angular  resolution,  which 
allows  us  to  use  it  for  the  search  of  inner  structure  of  a  light  current  within  a  beam  possessing 
optical  vortex. 

The  schematic  in  Fig.  1  shows  a  leaky-wave  coupler  using  two  prism  substrates  and  am 
immersion  liquid  between  them.  The  necessary  condition  for  the  liquid  refractive  index  to  be  less 
than  the  substrate  index  is  easily  realized.  Also,  it  makes  no  experimental  difficulties  to  achieve 
parallel  surfaces  of  the  film  with  a  thickness  of  several  wavelengths.  The  operation  in  the 
reflection-beam  detection  can  resolve  a  narrow  range  of  incident  angles  within  an  optical  beam 
with  diffraction-limited  divergence,  as  the  following  analysis  shows. 


Figure  1.  Sketch  of  the  waves  in  a  leaky-waveguide  device:  E0  is  the  incident  plane  wave,  ER  is 
the  reflected  wave,  ET  is  the  transmitted  wave.  Two  waves  Ex  and  E2  exist  within  the  liquid  film 
of  the  thickness  h  and  refractive  index  n2  between  the  prisms  with  refractive  index  nx.  X-axis  is 
parallel  to  the  film,  7-axis  is  perpendicular  to  the  plane  of  incidence,  Z-axis  is  directed  inside  the 


film. 


For  TE  polarized  plane  wave  with  amplitude  E0  incident  to  a  waveguide  liquid  film  with 
refractive  index  n2 ,  located  between  two  identical  substrates  of  refractive  index  nx,  the  boundary 
conditions  on  the  first  surface  of  a  film  (z  =  0)  will  look  as  follows: 


E0  +ER  -  E\  +E2’ 


n ,  cos  a 


9 


(i) 


where  E\  and  E2  are  the  amplitudes  of  waves  within  the  film,  ER  is  the  amplitude  of  the  reflected 
wave,  0Ci  and  oc2  are  the  angles  of  incidence  and  refraction,  as  indicated  in  Fig.  1. 

On  the  other  surface  (z  =  h)  the  boundary  conditions  are 


Ex  exp  \^hn2  cos  a2  0  E2  exp  b 


khn2  cos  a2 


E]  exp  \$hn2  cos a2  Qe2  exp  I Qlkhn2  cos a2  GJ 


GJ  Et  exp  \$hnx  cos  ax  Q 

-  Et  exp  I ^hnx  cos  ax  Q 


nx  cos  a. 


n2  cos  a2 


where  k  is  the  vavenumber  (in  vacuum),  E T  is  the  amplitude  of  the  transmitted  wave.  The  Snell 
low  relates  the  angles  0C\  and  a2\  ftisinoci  =  n2 sina2. 

Calculation  of  the  reflection  coefficient  R  =  |£r/£o|2  leads  to  the  result 


D  sin  (£/m2  cosa2) 

(3) 

•  2/77  \  n 2  5 

sin  (khn2  cosa2)  +  4p 

where 

„  n,n2  cos  a,  cos  a. 

'(4) 

B  = 

E  2  2 

Evidently,  the  transmission  coefficient  T=  1-7?.  The  smaller  j8,  the  sharper  is  the  resonance  peak 
in  transmission  and  the  gap  in  reflection,  which  correspond  to  the  leaky  mode  excitation.  The 
half- width  of  the  gap  can  be  estimated  to  be  Aax  =  2/J. 

The  resulting  reflection  coefficient  (intensity)  is  shown  in  Fig.  2  as  function  of  the  angle 
of  incidence  ax.  It  is  important  that  in  the  reflected  wave  the  condition  of  a  leaky  mode  excitation 
corresponds  to  a  zero  intensity.  Therefore,  for  an  ordinary  divergent  beam  (Gaussion  laser  mode, 
for  example)  a  narrow  dark  line  in  the  reflected  beam  cross-section,  oriented  perpendicularly  to 
the  plane  of  incidence,  will  appear  and  will  be  observable  in  the  near  field  as  well  as  in  the  far 
field.  Diffraction  spread-out  is  not  able  to  wash  out  the  zero-intensity  line  across  the  beam  due  to 
the  difference  in  phase  on  n  occurring  between  the  beam  sides  divided  by  the  dark  line.  This 
feature  is  very  important,  as  it  gives  a  sharp  notch  on  a  beam,  irrespectively  to  the  width  of  the 
gap  in  the  reflection  function. 

The  dark  line  in  the  reflected  beam  visualize  a  narrow  range  of  angles  of  incidence  within 
a  beam,  corresponding  to  the  effective  leaky  mode  excitation.  Having  in  hands  this  sensitive 
instrument,  we  can  apply  it  for  the  analysis  of  inner  structure  of  a  helicoidal-wavefront  optical 
beam,  or  a  beam  possessing  axial  optical  vortex  [2]. 


Angle  of  incidence  (deg) 

Fig.  2.  Calculated  angular  dependence  of  the  reflection  from  a  leaky  waveguide 
(glass/glicerine/glass).  The  angle  70.5°  corresponds  to  the  total  internal  reflection. 


To  analize  the  shape  of  helicoidal  wave  front  inside  the  angle  of  diffraction  limit  the 
“ray”  model  of  OV  beam  was  used.  Figure  3  represents  a  calculation  of  the  line  being  a  cross- 
section  of  the  ray  family  which  have  in  YZ  plane  zero  inclination  with  respect  to  the  beam  axis  (8 
=  0).  All  these  rays  will  excite  a  waveguide  mode,  what  can  be  detected  as  in  the  transmitted 
beam,  as  in  the  reflected  beam.  Observation  of  the  dark  line  in  the  reflected  beam  will  give 
therefore  a  visual  picture  of  the  loci  of  rays  with  equal  angle  of  incidence  in  the  YZ  plane. 


(a)  (b) 


Fig.  3.  Calculation  of  the  cross-section  of  the  loci  of  ray  family  with  zero  inclination  with  respect 
to  the  beam  axis,  in  the  plane  of  incidence.  In  calculations,  we  let  k  =  10  pm'1,  R  =  1.5T06  pm, 
the  radial  scale  is  given  in  micrometers  and  the  azimuth  scale  in  degrees.  For  zero  detuning  (8  = 
0)  all  these  rays  will  excite  a  propagating  leaky  mode  in  the  waveguide,  thus  producing  a  dark 
line  in  the  reflected  beam  (a).  (b)The  shape  of  the  dark  line  in  the  reflected  beam  for  constant  ray 
inclination  with  respect  to  the  beam  axis,  which  compensates  the  beam  angular  detuning  8  =  - 
4-1  O'4. 


The  shape  of  the  line  strongly  depends  on  the  wave-front  curvature  and  the  detuning  5. 
Non-zero  detuning  shifts  the  dark  line  from  the  beam  axis.  The  smaller  is  the  spherical  component 
of  the  curvature,  the  sharper  is  the  influence  of  the  helical  componentgives  an  For  example,  at  R  = 
°o?  and  8  =  5T0'5:  the  line  turns  now  to  a  closed  loop.  This  situation  corresponds  to  the  focused 
laser  beam,  when  the  incident  beam  enters  the  waveguide  just  in  the  waist.  For  another  sign  of  the 
detuning,  the  loop  will  appear  in  the  upper  half  plane.  Zero  detuning  in  this  case  corresponds  to 
the  horizontal  line  through  the  beam  centre. 

Two  identical  60°  glass  prisms  with  refractive  index  n\  =  1.515  were  used  to  compose  a 
leaky-mode  waveguide.  A  drop  of  glycerine  (n2  =  1.43)  placed  between  prisms  produced  a 
parallel  film  after  the  stressing  the  prism  faces.  The  thickness  of  the  layer  was  estimated  to  be 
about 

The  quality  of  just  obtained  waveguide  was  tested  with  an  expanded  (divergent)  laser 
beam  (Gaussian  mode)  from  He-Ne  laser.  A  set  of  vertical  narrow  dark  lines  appears  in  the 
reflected  beam.  With  non-expanded  beam  from  the  laser  output,  the  resonance  conditon  of  leaky 
mode  excitation  can  be  detected  by  observation  of  one  dark  line  crossing  the  beam  (Fig.  4a). 


(b) 


Fig.  4.  Dark  line  crossing  the  Gaussian  beam  reflected  from  the  leaky  waveguide,  zero  angular 
detuning  (a),  the  same  for  an  incident  OV  beam,  non-zero  detuning  (b),  zero  detuning  (c). 


To  create  an  OV  beam,  we  used  a  usual  technique  of  laser  beam  (TEM  mode  from  HeNe 
laser)  diffraction  by  a  computer-synthesized  hologram  [2]. 

The  divergence  of  the  OV  beam  obtained  in  the  first  order  of  diffraction  (m  =  1)  was 
measured  to  be  about  2°.  The  OV  beam  was  directed  to  the  prism-film  coupler,  which  was 
mounted  in  a  table  rotating  with  high  precision. 

In  order  to  verify  in  one  experimental  scheme  the  predictions  of  the  theoretical  analysis 
both  for  convergent  and  divergent  beams,  we  used  a  lens  with  focal  length  30  cm  and  gradually 
move  the  prism  coupler  before  and  behind  the  focal  plane. 

The  images  of  the  reflected  beam  taken  by  CCD  camera  are  shown  in  Fig.  4b, c.  The  angle 
of  incidence  was  chosen  to  observe  the  first  mode  resonance  in  the  reflected  beam.  For  non-zero 
angular  detuning,  there  is  a  slightly  bent  dark  line  crossing  the  beam,  in  accordance  with  the 
geometrical  optics  analysis  (Fig.  4b). 

There  is  noticeable  deviation  from  the  expectation  folloing  from  the  geometrical  optics 
description  in  the  case  of  zero  detuning.  Just  in  the  vicinity  of  the  OV  axis  we  have  detected  a  tear 
of  the  dark  line  and  disappearence  of  the  central  zero-intensity  point  which  is  the  OV  core  (Fig. 
4c),  instead  of  the  calculated  continuous  dark  line  crossing  the  OV  core  (Fig. 3). 

As  a  summary,  the  used  leaky-mode  waveguide  was  shown  to  be  very  sensitive 
instrument  for  a  detail  investigation  of  the  structure  of  an  OV  beam  with  diffraction-limited 
divergence.  The  reason  is  the  multiple-beam  interference  within  a  liquid  film,  producing  high 
angular  resolution,  especially  for  low-order  modes  (the  quality  of  the  wavegide  is  higher  for  them 
due  to  the  higher  reflection  on  the  boundaries). 
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2.8  Second  harmonic  generation  of  OV  beams  with  fractional 

topological  charge 


We  have  shown  that  there  is  a  possibility  to  create  optical  beams  which  possess 
optical  vortices  with  fractional  topological  charge  by  a  holographic  technique  [1]. 

Now  we  report  an  elaborated  technique  and  the  experimental  results.  We  also 
discuss  some  nontrivial  effects  of  second-harmonic  generation  (SHG)  of  OV  beams 
with  fractional  topological  charges. 

One  of  the  most  important  properties  of  an  optical  vortex  is  the  circulation  of  phase 
around  its  axis,  in  the  form  [2] 

0(z  =  const )  =  777  (p ,  (1) 

where  (p  is  the  azimuth  angle  and  m  is  the  topological  charge  of  an  OV.  This  transversal  phase 
structure  is  responsible  for  the  specific  interference  pattern  produced  with  an  ordinary  tilted 
reference  wave  (Fig. la).  Instead  of  a  usual  periodic  fringes,  the  splitting  of  one  fringe  appears 
into  m  new  fringes  [3].  The  interference  with  a  coaxial  plane  reference  wave  produces  spiral 
fringes,  as  shown  in  Fig.  lb. 


Fig.  1.  Pattern  of  interference  between  charge  1  OV  beam  and  tilted  plane  reference  wave  (a)  and 
between  charge  4  OV  beam  and  coaxial  plane  reference  wave  (b). 

The  idea  of  OV  beam  formation  by  use  of  diffraction  of  an  ordinary  wave  on  a  computer¬ 
generated  grating  is  based  on  the  holographic  principle:  a  readout  beam  restores  the  wave,  which 
participated  in  the  hologram  recording.  Instead  of  writing  a  hologram  with  two  actual  optical 
waves,  it  is  enough  to  calculate  the  interference  pattern  numerically  and  print  the  picture  in  black 
and  white  or  gray  scale.  Then  the  picture  after  reduction  of  the  transverse  dimensions  serves  as  an 
amplitude  grating  producing  OV  beams  in  diffraction  orders  [4].  Examples  of  binary  gratings  are 
shown  in  Fig.  2. 


(a)  (b) 


Fig.  2.  The  computer-generated  patterns  for  gratings  which  are  able  to  produce  in  the  first 
diffraction  order  OV  beams  with  topological  charges  m  =1  (a)  and  m  =  5  (b). 


The  holographic  principle  can  be  applied  also  for  the  production  of  more  complicated  beams  with 
specific  features.  For  the  task  of  OV  beams  with  fractional  topological  charge  creation,  it  is 
necessary  to  calculate  the  interference  pattern  which  corresponds  to  the  phase  dependence, 
according  to  Eq.(l),  with  the  fractional  charge  m.  An  example  of  a  binary  grating  for  m  =  lA  is 
shown  in  Fig.  3a  [1].  A  half-cut  of  the  fringes  appears  along  a  line  (p  =  0,  due  to  the  phase  step  on 
7t  on  this  line.  However,  this  method  is  not  unique.  We  propose  now  more  interesting  technique, 
which  is  based  on  the  production  of  subharmonic  diffraction  orders.  Let  us  suppose  that  the 
grating  function  has  the  form 

T(p,(p)  =  T0  cos (p + <p  +  (p0 GJa2  cosj^COS<^  +  <^>  +  <^0 

With  the  condition 

2 

a 7 

ao=ai+  —  >  (3) 

8  ax 

formula  (2)  gives  minimum  transmission  Tmin  =  0,  and  periodically  alterating  maximum 
transmission  between  Tmaxl  =  7b(tf0  +  ^ i  +^2)  and  Tmax2  =  To(a0  +  a\-  a2).  The  structure  of  the 


transmission  function  is  shown  in  Fig.  3.  The  choice  of  <p0  =  37E/2  permits  to  obtain  a  symmetric 
function.  Figure  3  shows  the  transmission  function. 


Fig.  3.  Mesh  of  the  transmission  function  of  a  grating  with  m  =  lA  subharmonic,  according  to 
Eq.(2). 

The  diffraction  by  the  grating  with  modulated  contrast  will  produce,  except  ±1  diffraction 
orders,  which  corresponds  to  the  fundamental  period,  a  pair  of  subharmonic  orders,  which 
will  be  spatially  located  on  the  bisectrix  between  the  ±1  diffraction  orders.  In  a  similar  way 
we  can  construct  gratings  which  can  produce  any  necessary  fractional  charge  of  an  OV. 

The  rectangular-profile  pattern  for  a  computer-generated  diffraction  grating  was  calculated  on 
the  basis  of  the  transmission  function 
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integer  n  =  1/m,  A  is  the  grating  spacing. 


The  example  for  a  grating  pattern  with  n  =  2  (m  =  V2)  is  shown  in  Fig.4c. 


Fig.4.  Three  examples  of  computer-generated  gratings,  suitable  for  the  production  of  OV  beams 
with  fractional  topological  charge  m  =  V2.  (a)  Binary  grating  with  a  half-cut  of  fringes  along  a 
horizontal  line;  (b)  Binary  grating  with  a  half-cut  along  a  vertical  line;  (c)  Grating  with 
periodically  modulated  contrast 

The  proposed  techniqe  is  easily  applied  for  the  creation  of  any  subharmonic  beams  with 
fractional  charges  m  =  1  In,  2 In,  Figure  5  represents  the  diffraction  orders  for  the 

computer-synthesized  grating  shown  in  Fig.  4b.  A  half-cut  oriented  roughly  vertically  appears  in 
the  subharmonic  beam  just  behind  the  grating,  and  then  gradually  transforms  to  a  single  OV  in  the 
far  field.  Subharmonic  beam  with  the  topological  charge  3/2  has  more  complicated  structure  in 
the  near  field,  and  possesses  two  isolated  OVs  in  the  far  field. 


Fig.  5.  Diffraction  orders  in  the  near  field,  from  left  to  the  right:  zero,  subharmonic  1/2,  first, 
subharmonic  3/2. 

Figure  6  shows  the  diffraction  orders  for  the  grating  produced  to  create  1/3  subharmonic, 
as  well  as  2/3  subharmonic  beams.  In  the  far  field  both  subharmonics  (1/3  and  2/3)  have  one 
single-charged  OV. 


Fig.  6.  Diffraction  orders,  from  left  to  right:  zero,  subharmonic  1/3,  subharmonic  2/3,  first. 

Both  gratings  (shown  in  Figs  5  and  6)  give  a  single-charged  axial  OV  in  the  first 
fundamental  order,  which  remains  the  same  up  to  the  far  field. 

As  it  was  established  earlier,  OV  beams  with  fractional  topological  charges  are  not  stable 
along  the  propagation  path,  in  contrast  to  an  axial  charge  1  vortex  (e.g.  Laguerre-Gaussian  mode 
LG1  o)-  A  decomposition  of  the  diffracted  beam  to  the  set  of  LG  modes  shows  an  immidiate  decay 
of  the  initial  half-cut  of  the  wave  front  to  an  infinite  set  of  individual  vortices,  which  gradually 
leabe  the  beam  along  its  propagation  path.  In  the  far  field,  the  number  of  vortices  is  equal  to 
mod(m). 


(C)  (d) 


Fig. 7.  Evolution  of  an  OV  beam  with  fractional  charge  1/2  in  the  near  field,  (a)  The  Vi 
subharmonic  diffraction  order  on  a  distance  16  cm  from  the  grating  and  (b)  the  corresponding 
interferogram.  The  vertical  dark  line  crosses  the  beam,  there  is  a  phase  step  seen  in  the 
interferogram.  Individual  vortices  are  not  resolved,  (c)  The  diffracted  beam  on  the  distance  100 
cm  from  the  grating.  The  central  isolated  vortex  is  seen,  several  peripheral  vortices  can  be  also 
found  in  the  interference  pattern  (d). 

With  respect  to  the  process  of  second-harmonic  generation  (SHG)  with  OV  beams,  an 
interesting  optical  paradox  appears.  The  phase  doubling,  which  occurs  in  the  SHG  and  leads  to 
the  doubling  of  the  topological  charge  of  the  OV  [1,5],  will  create  a  single-charged  OV  from  the 
fractional  charge  Vi  in  the  fundamental  frequency  (FF)  beam.  From  the  other  hand,  if  the 
nonlinear  frequency-doubling  element  will  be  located  far  enough  from  the  grating,  the  incident 
beam  will  contain  a  single  off-axis  OV,  and  two  off-axis  OVs  are  expected  to  appear  in  the  SH 
beam  [1,5].  In  the  intermediate  situation,  a  set  of  vortices  existing  in  the  FF  beam  will  produce  a 
doubled  set  of  OVs  in  the  SH  beam.  The  evolution  of  this  set  of  OVs  will  take  place  in  the  SH 
beam,  and  the  result  is  not  evident  in  advance. 
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Chapter  3 

Multichannel  Optical  Manipulator 
Based  on  OV  Beam 

3.1  Introduction 

The  physical  background  of  manipulations  (translation,  rotation,  arrangement,  etc.) 
by  OVs  with  captured  micron-size  particles  is  considered.  Manipulation  of  biological  objects 
allows  to  carry  out  biological  reactions  between  selected  components  or  to  accelerate  natural 
reactions  by  means  of  closing  reacting  components. 

In  1974  Ashkin  and  Dziedzic  [1]  showed  that  TEM*0i  “doughnut”  laser  mode  could 
be  used  to  levitate  various  types  of  spheres.  They  described  levitation  of  hollow  transparent 
and  metallic  spheres  and  also  solid  transparent  spheres.  Solid  transparent  spheres  levitated 
within  a  single  vertically  upward  propagating  doughnut  beam  while  hollow  transparent  and 
metallic  spheres  was  levitating  under  action  of  a  combination  of  two  counter-propagating 
beams  and  one  vertically  upward  beam.  Physical  base  of  optical  trapping  arises  from 
Lebedev’s  work  on  measurement  of  radiation  pressure  in  1901.  Firstly,  a  particle  exchanges 
the  momentum  with  the  light  field.  Next,  the  local  flow  of  momentum  is  configured  in  such  a 
way  so  that  an  equilibrium  point  exists.  The  mechanism  of  particle  trapping  with  single  OV 
was  studied  by  Heckenberg  and  collaborators  in  detail  [2]. 

3.2  Computer  simulation  of  the  multichannel  optical  manipulator 

Contrasting  to  the  aforeciting  researches  on  controlled  3D  and  2D  trapping,  we 
propose  the  theoretical  base  for  trapping  of  two  microparticles  (or  group  of  them)  with  a 
possibility  of  simultaneous  manipulation  of  them.  All  manipulations  are  possible  in  each 
cross-section  plane  of  the  beam. 

We  constructed  the  beam  with  two  embedded  OVs  of  opposite  signs  (OV  “dipole” 
OVD)  disposed  on  some  distance  from  the  beam  centre  [3].  Such  beam  can  be  produced  easy 
by  superposing  of  several  HG  modes  or  the  combination  of  the  HG  and  the  LG  modes  with  a 
necessary  phase  shift.  The  following  mode  composition  can  be  used: 

LGo1  +  2i  HGo1  +  HG0°  e'¥  «  ‘Vortices  Dipole’ ,  ( 1 ) 

where  \| /  is  the  phase  parameter,  introduced  to  control  the  OVD  size.  The  distribution  of  the 
intensity  in  the  cross-section  of  OVD  is  depicted  in  Fig.  1.  Each  dark  spot  is  the  OV  core. 
Microparticles  within  the  OV  cores  happen  to  be  captured  radially.  Thus,  we  have  a  two- 
channel  trap.  Each  OV  in  OVD  has  the  local  intensity  and  phase  structure  near  the  OV  axis 
similar  to  a  single  OV. 


(a)  (b)  (c) 


Fig.  1.  Intensity  distribution  in  the  cross-section  of  the  Gaussian  beam  with  embedded  OVD: 
(a)  the  initial  state,  (b)  close  to  the  annihilation  state  and  (c)  after  the  annihilation. 


The  ability  of  OVD  to  change  the  size  and  finally  to  annihilate  under  incursion  of 
some  phase  shift  of  the  HG0o  mode  is  the  base  of  the  multichannel  OV  manipulator.  We  use 
the  phase  shift  \| /  to  control  the  OVD  size.  Figure  2  presents  the  OVD  size  vs  phase  shift.  In 
experiment  the  phase  variation  can  be  easily  controlled. 


Fig.  2.  OVD  size  vs  phase  shift  \| /:  (a)  one  Rayleigh  distance  from  the  waist;  (b)  5  Rayleigh 
distances,  (c)  10  Rayleigh  distances,  (d)  20  Rayleigh  distances.  The  phase  scale  (horizontal 
axis)  is  measured  in  radians,  multiplied  by  100.  Vertical  axis:  the  OVD  size  is  presented  in 
units  normalized  on  the  current  transversal  size  of  the  Gaussian  beam. 


Next  our  step  is  the  tree-channel  manipulator.  It  bases  on  the  following  principle, 
with  only  two  LG  modes  used  to  compose  the  manipulator  (all  three  vortices  have  the  same 
sign): 

LG03+A  LG0°  exp(i\|/)  «  Three-channel  manipulator’,  (2) 


where  A  and  \| /  are  the  control  parameters.  As  for  the  case  of  two-channel  manipulator,  we 
introduce  here  the  phase  \|/  in  Gaussian  LG0°  mode.  Besides,  a  new  control  parameter  A 
appears,  which  is  responsible  for  the  distancing  of  all  channels  radially.  The  larger  is  A,  the 
larger  is  the  distance  of  the  OV  from  the  centre  of  the  beam.  When  A  equals  to  zero,  the 
channels  are  not  separated. 

Unlike  previous  case  of  two-channel  OVD  system,  \|/  is  the  azimuth  parameter 
responsible  for  the  rotation  of  all  channels.  The  increasing  \|/  to  2n  results  in  the  one  full  turn 
of  the  whole  system  of  channels.  Decreasing  of  \|/  leads  to  the  pattern  rotation  in  opposite 
direction. 


(a) 


(b) 


Fig. 3.  (a)  radial  motion  of  the  channels  while  A  is  increasing  and  circular  motion  when  \| /  is 
increasing;  (b)  spiral  channels  motion  while  A  and  \|/  are  increasing  simultaneously. 

3.3  Experimental  realization:  analysis  of  possibilities  of  the 
manipulator 

The  operating  parameters  were  chosen  widely  by  comfortable  way  of  realisation  in 
experiment.  The  variation  of  phase  \|/  can  be  produced  with  a  mirror  on  the  piezo-element  put 
in  optical  setup.  The  use  of  piezo-elements  for  small  space  displacement  is  widely  used 
technology  in  optics. 

The  parameter  being  the  LG0°  mode  amplitude  can  be  changed  by  many  ways.  The 
simplest  one  is  the  introducing  of  an  attenuator  on  the  way  of  the  beam.  In  our  opinion,  it  is 
better  to  use  the  dielectric  mirror,  which  possesses  angular  dependence  of  reflectivity.  By  the 
tuning  of  such  mirror  we  are  able  to  operate  the  intensity  of  the  passing  beam. 

The  principle  used  for  the  constructing  of  three-channel  manipulator  can  be  extended 
to  the  construction  of  a  manipulator  with  any  number  of  channels.  Two  modes  are  enough  to 
construct  these  manipulators:  LG0m  and  LG0°  modes,  where  \m\  is  number  of  the  channels.  The 
advantage  of  such  approach  is  that  the  number  of  operating  parameters  for  manipulation 
remains  two. 
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Chapter  4 

Optical- vortex  fiber  sensors  (OVFS)  with  cw  and  pulse  operation 


Summary 

The  theoretical  and  experimental  study  of  optical-vortex  fiber  sensors  (OVFS)  and  their 
assessment  are  performed  The  computer  simulation  of  applicability  of  short-pulse  lasers  in  OVFS 
is  carried  out. 

As  a  basis  of  OVFS,  the  physical  approach  involving  the  analysis  of  a  light  beam 
propagation  along  a  OV-generated  unit,  OV-guided  unit  and  OV-detection  unit  is  developed.  A 
computer  processing  technique  of  an  output  interference  image  was  elaborated. 

OV  decay  in  a  fiber  was  considered.  Exploitation  of  optical  fibers  preserving  a  guided 
OV  in  supersensitive  fiber-optical  devices  for  a  monitoring  both  temperature  parameters  and 
inertial  navigation  systems  is  discussed. 

4.1  Analytical  review  and  statement  of  the  problem 

The  special  properties  of  OVs  enable  us  to  develop  new  optical  measurement  and 
research  methods  [1-4].  In  particular,  the  studies  on  the  OV  phenomena  can  also  extend  our 
understanding  of  the  topological  manifestations  in  different  fiber-optical  effects  and,  in  turn,  new 
fiber-optical  devices. 

Interferometric  fiber-optical  sensors  have  been  the  subject  of  considerations  for  last 
twenty  years  [5].  The  application  of  OVs  as  probe  and  signal  waves  may  stimulate  new  success  in 
this  area. 

Although  some  types  of  fiber-optical  sensors  on  the  base  of  high-order  fiber  modes  have 
been  described  before  [6-8,12-14],  their  use  for  OV-based  sensor  are  presented  here  for  the  first 
time. 

Consider  at  first  the  main  features  of  two-  and  one-arm  single-mode  fiber  interferometers. 

4.1.1  Two-arms  fiber  interferometer 


Let  the  length  of  the  reference  arm  be  Li  and  the  signal  arm  (fiber)  L2.  The  reference  arm 
is  thermally  stabilized,  and  the  signal  arm  is  affected  by  the  temperature  variations.  Hence, 
relative  phase  changes  are  restricted  to  the  variations  of  the  signal  arm  parameters.  Therefore  the 
phase  in  the  reference  arm  is  constant,  and  in  the  signal  arm  is 
2 n 

<Z>2=  —  nL2,  (1) 

A 

where  fT  is  the  effective  refractive  index  of  the  fiber  for  the  fundamental  H  En  mode  of  the  fiber, 

X  is  the  vacuum  wavelength  of  the  used  source. 

The  temperature  sensitivity  of  the  system  arises  from  thermal  expansion  and  from  the 
temperature  dependence  of  the  refractive  index  in  the  signal  arm: 


1  d<&2  2k  f  n  d  L2  dn  ^ 

17  dr  “X  l  ar  ar  ’ 

2  v^2 

where  T  is  the  absolute  temperature.  Typically  [7] 


(2) 

-^TocSxlO-'/r1,  ^-oclO-'/r1,  hence 


the  temperature  dependence  of  the  refractive  index  dominates  and 


1  d<t>  ,  „rjad 
- oc  100— 

LdT  Km 


for  l  =0.63 


pm. 

A  number  of  data  processing  schemes  for  the  fiber  interferometer  have  been  developed 
[5,7]  which  yield  phase  resolution  in  the  lfirdd  +  lmrad  range.  The  particular  value  of  the 


temperature  sensitivity  depends  on  an  optical  fiber  structure.  The  birefringence  single-mode  fibers 
are  affected  by  the  most  temperature  action  on  account  of  strong  internal  tensions.  For  example, 
the  birefringence  fiber  with  the  beat  length  up  to  1  mm  gives  the  temperature  sensitivity  more 

than  200  [7]. 

Km 

4.1.2  One-arm  interferometer 

Other  type  fiber  optical  interferometers  are  one-arm  devices.  In  contrast  to  the  two-arm 
interferometers  those  are  extremely  simple  constructions.  It  requires  only  two  polarizers  for  the 
splitting  and  reconstruction  of  the  two  interfering  modes.  A  particularly  advantageous  sensor  in 
this  respect  is  a  fiber-optical  interferometer  employing  a  fiber  with  high  internal  lateral  stresses 
[9].  This  fiber  is  highly  birefringent,  and  two  eigen-modes  are  the  fast  and  slow  H  En  modes  of 
mutually  orthogonal  linear  polarizations.  The  temperature  sensitivity  is  due  to  the  thermal 
expansion,  which  affects  the  magnitude  of  the  internal  lateral  stress  and,  therefore,  the  fiber 
birefringence.  Simultaneously,  this  fiber  is  fairly  insensitive  to  isotropic  pressure  variations. 

In  all  interferometric  sensors  discussed  here,  the  phase  response  <3>  is  proportional  to  the 
length  L  of  that  section  of  a  fiber,  which  is  exposed  to  the  external  influence  Fto  be  sensed.  The 
influence  F  may  be  caused  by  a  variation  either  of  temperature  T  or  of  other  ambient 
thermodynamics  values.  The  phase  delay  of  modes  in  a  fiber  of  length  L  is 

AO  =  k  (fij  —  /T2 )  L  —  kAn  L  ,  where  AfT  is  the  difference  of  the  mode  effective  indices.  The 
phase  response  becomes 

0) 

L  dF  dF  LdF 

This  is  a  /--sensitivity  of  the  one-arm  interferometer.  If  R  is  a  birefringence,  we  have  from  Eq.  (3) 

d(Af))  d(AfT)  d(A  fT)dR 

dF  dF  dR  dF  J 


The  temperature  sensitivity  for  a  stress-induced  single-mode  fiber  may  be  expressed  as 
[9]: 

1  d  (AO)  k  AfTdfT  dL  A  nir ) 

- .  L  +  AfT  ,+L 

L  n  dT  dT  T-L 


L  dT 


(5) 


J 


where  Ts  is  the  softening  temperature  of  the  inner  cladding.  So  the  temperature  sensitivity  for  the 
birefringent  single-mode  fiber  with  the  beat  length  A  =  3. 2  /DfD  or  AfT  —  2  X  10-4  is 
1  d  (AO)  rod 

—  2.48 - ,  while  the  temperature  sensitivity  of  the  two-arm  interferometer  with  the 


L  dT 
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same  fiber  [9]  becomes 


1  d(  AO) 
T  dT 


- 124 


rad 

Km 


In  order  to  increase  the  temperature  sensitivity  of  the  one-arm  interferometers  one  makes  use 
of  different  special  adjustments  for  a  fiber.  To  the  end,  for  example,  an  isotropic  single-mode 
fiber  is  coiled  under  some  tension  around  a  hollow  cylinder  having  a  high  thermal  expansion 
coefficient  [10],  In  such  a  case  the  fiber  becomes  equivalent  to  a  retardation  plate  with 
birefringence  depending  on  the  outer  radius.  When  the  temperature  is  increased  the  two  principal 
indices  of  the  fiber  change  because  of  thermooptic  effect  and  elastooptic  effect  caused  by  the 
thermal  expansion  of  the  cylinder.  As  a  consequence  the  phase  shift  of  the  eigen-modes  are 


2 


has 


modified  by  different  amounts.  In  this  way  the  temperature  sensitivity 


1  d  (A<E>)  rad 
L  dJ  I  U Km 


been  reached. 

The  other  way  of  the  temperature  sensitivity  increase  is  concerned  with  the  spectral 
registration  of  a  fiber  output  signal  [1 1],  If  a  white  light  excites  the  fiber  the  modes  with  different 
wavelengths  has  a  different  temperature  response.  This  high  accuracy  method  for  spectral  data 

processing  allows  achieving  a  temperature  resolution  of  3x  10  3  K  for  a  wide  temperature  range. 


4.1.3  Dual-mode  fiber  sensor 

In  the  works  [6,12-14]  the  dual-mode  acoustic  and  temperature  sensors  on  the  base  of  a 
multimode  fiber  are  discussed.  The  multimode  fiber  turns  out  to  enable  to  increase  the  sensitivity 
of  the  fiber-optical  sensors.  Two  modes  of  the  fiber  ( LP01  and  LPn  )  form  two  signal  arms  and 
produce  differential  phase  change  between  the  beams  propagating  in  the  two  modes. 


(a)  (b) 

Fig.l.  Intensity  distribution  of  LP0l  (a)  and  L  Pn  (b)  modes  at  the  fiber  cross-section  [6]. 

The  fiber  used  in  the  experiments  of  the  work  [6]  was  a  conventional  quasi-step  index 
communications  type  fiber  designed  to  be  single-mode  at  1.3  |Lim.  The  experiments  were 
performed  in  the  1  |Lim  wavelength  region  where  the  fiber  guides  two  modes,  according  to  [6]. 
The  polarization  state  of  the  two  modes  was  the  same.  As  it  can  be  seen  from  Fig.  lb,  the  intensity 
structure  of  the  LPn  mode  represents  two  lobes  because  of  the  71  -  phase  difference  between 
them.  The  phase  variations  take  place  due  to  the  different  phase  velocities  of  LP01  and  LPn 
modes.  The  variation  in  the  phase  difference  can  be  seen  as  a  transfer  of  power  between  two 
lobes.  The  simplest  detection  scheme  selects,  spatially,  a  single  lobe  from  the  fiber  near  field,  and 
as  the  phase  difference  is  varied,  the  detected  intensity  oscillates  through  a  series  of  maxima  and 

minima.  This  method  of  the  signal  registration  gives  only  0.18-^^,  while  the  spectral 
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rad 

registration  method  raises  the  temperature  sensitivity  up  to  150  -  [6]  (the  sensor  with  the 

Km 

spectral  registration  method  is  considered  Section  4.3.3). 

4.1.4  Statement  of  the  problem 


3 


Fiber-optical  sensors  have  been  a  research  field  in  their  own  right  since  1975.  Some  sensors 
have  been  the  subject  of  intense  research,  particularly  for  hydrophone  and  gyroscope  applications, 
and  some  extremely  complex  and  highly  sensitive  devices  have  been  engineered. 

Nevertheless,  there  are  a  large  areas  which  have  received  relatively  little  attention  and  there 
is  much  scope  for  the  development  of  novel  and  improved  sensing  devices,  in  particular,  the 
fiber-optical  sensing  on  the  base  of  multimode  fibers. 

Indeed,  development  of  those  highly  sensitive  schemes  as  two-arms  interferometers  (for 
example,  Mach-Zehnder  interferometer  configuration)  is  faced  with  great  difficulties.  The  main 
problems  with  this  technique  are  connected  with  occasional  phase  shift  caused  by  physical 
variations  other  than  ones  to  be  measured.  At  the  same  time,  the  two-arms  interferometric  sensors 

rad 

of  temperature,  for  example,  have  the  temperature  sensibility  more  than  100 - . 

Km 

These  problems  are  partially  solved  by  means  of  utilizing  of  the  one-arm  interferometers, 
in  particular,  the  devices  on  the  base  of  multimode  fibers.  These  temperature  sensors  are  able  to 

rad 

raise  the  temperature  sensitivity  up  to  200 -  and  more  if  the  spectral  registration  of  signals  is 

Km 

used.  However  the  intensity  registration  conceals  some  problems.  The  main  one  is  concerned  with 
noise  sources  such  as  laser  instability  or  laser  intensity  fluctuations. 

Moreover,  the  use  of  a  multimode  fiber  assumes  a  monomode  or  duel-mode  regime,  i.e. 
one  or  two  modes  must  propagate  along  the  fiber  only.  In  turn,  this  tough  requirement  is  closely 
associated  both  with  the  excitation  conditions  of  the  fiber  and  concordance  of  mode  polarizations. 
So,  authors  of  the  work  [6]  claim  that  the  modes  L  PQ1  and  L  Pn  transmit  through  the  fiber  sensor 
as  the  eigen-modes,  and  these  modes  have  the  same  polarization  state.  However,  the  LPn  mode 
is  known  [15]  to  be  no  the  fiber  eigen-mode  but  to  be  either  the  sum  of  TE01  and  H  or 

TM  01  and  HE  fiber  eigen-modes.  Moreover,  their  polarization  state  changes  along  the  fiber 
and  the  linear  polarization  is  recovered  at  certain  planes  of  the  fiber  cross-section  only.  Hence,  the 
dual-mode  fiber  sensor  is  an  instable  device  with  a  strong  mode  noise:  any  slight  fiber 
perturbations  cause  the  experimental  mistakes. 

In  order  to  avoid  or  to  reduce  the  intensity  instability  it  is  necessary  to  interchange  the 
registration  method  in  principle  by  another  one.  The  mode  noise  may  be  reduced  by  means  of  the 
excitation  of  a  relatively  more  stable  mode  superposition  such  as  circularly  polarized  fields. 

The  OV  fiber  modes  have  these  properties.  On  the  one  hand,  the  singular  beams  are  the 
circularly  polarized  eigen-modes,  which  propagate  along  the  circular  optical  fiber  without 
changes  of  their  field  structure  [16,17].  On  the  other  hand,  the  interference  pattern  of  the  optical 
vortex  and  fundamental  Gaussian  beam  represents  the  spiral-like  image  rotating  around  its  axis 
due  to  the  phase  difference  between  the  signal  and  reference  beams  [18,19].  Besides,  the  rotation 
angle  of  the  spiral  image  is  directly  proportional  to  the  phase  difference  between  the  beams.  In 
this  case  the  registration  process  becomes  different:  the  data  processing  is  replaced  by  the  image 
processing.  This  enables  us  to  avoid  the  influence  of  laser  intensity  fluctuations  and  to  raise  the 
sensitivity  of  the  fiber-optical  device. 

On  the  base  of  the  above  listed  conclusions  the  follow  problems  of  investigations  were 

studied: 

1)  Theoretical  and  experimental  study  of  OVFS. 

2)  Consideration  and  computer  simulation  of  applicability  of  short-pulse  lasers  in  OVFS. 
Theoretical  and  experimental  solution  of  these  problems  will  be  done  in  the  following  sections. 
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4.2  Optical-vortex  fiber  sensors 


The  review  performed  above  has  showed  numerous  advantages  of  OVFS  on  the  base  of 
low-mode  and  single  mode  fibers  as  against  other  technologies.  Such  their  properties  as  high 
sensitivity,  resistivity  to  external  electromagnetic  perturbations,  geometrical  versatility,  safety  in 
hazardous  environment  and  so  on  make  them  irreplaceable  devices  for  quite  a  number  of 
applications.  We  shall  consider  sensors  based  on  low-mode  fibers  and  further  restrict  our  attention 
to  those  worked  into  single- vortex  regime. 

4.2.1  Guided  optical  vortices  and  devices  on  their  base 
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OV  transmittion  through  inhomogeneous  medium  of  an  optical  fiber  has  essential 
differences  with  those  propagation  in  free  space.  Without  going  into  detail  discussed  in  the  works 
[1-6]  at  full  length,  we  outline  general  traits  of  these  guided  field  forms. 

At  once  it  is  worthy  to  note  that  guided  OVs  are  the  main  form  of  a  field  existing  in 
circular  optical  fibers.  All  other  eigen-modes  of  a  fiber  can  be  presented  as  an  OVs  superposition 
since  the  guided  OVs  by  themselves  are  the  eigen-modes  written  in  a  circularly  polarized  basis 
[1,2,5].  There  is  exception  of  this  rule:  azimuthally-symmetric  TE  and  TM  modes  can  not  be 
written  in  the  circularly  polarized  basis  because  of  their  specific  field  structure  [2,3].  The  fact  is 
that  the  guided  OV  is  essentially  vector  field  structure,  due  to  impossibility  to  separate  its 
topological  charge  Z  and  circular  polarization  state  characterized  by  a  helicity  G  =  ±1  (the  sign  of 
the  circular  polarization).  This  OV  property  is  a  consequence  of  a  spin-orbit  interaction  in  the  OV 
field  [4].  The  difference  between  OV  and  TE  and  TM  modes  manifests  itself  as  polarization 

corrections  8(3  to  a  scalar  part  of  the  propagation  constant  of  a  mode  (3  :  (3  =  (3  +  8(3  .  The 
polarization  corrections  for  TE  and  TM  modes  are  strongly  differed  so  that  their  superposition 
TM  ±iTE  is  unstable  [2].  Although  this  superposition  forms  vortex-like  structures  at  some  fiber 
lengths,  in  the  intermediate  distances,  there  is  a  mode  mixture  without  definite  topological  charge 
and  helicity  (the  polarization  state).  This  mode  superposition  is  called  sometimes  as  IV  mode 
(instable  vortex).  In  contrast  to  the  IV  mode,  the  guided  OV  called  CV  mode  (circularly  polarized 

vortex)  in  any  point  has  definite  state  with  both  positive  topological  charge  and  helicity  |+  l,+l) 

or  negative  ones:  |—  l,-l) .  There  are,  as  a  rule,  one  CV  mode  and  one  IV  mode  in  a  low-mode 
circular  fiber  (1=0,1).  In  an  anisotropic  fiber  such  as  a  strongly  elliptic-core  fiber  or  a  stress- 
induced  core  fiber,  the  CV  mode  in  a  pure  state  cannot  exist  [6].  Nevertheless,  a  slightly 
birefringent  fiber  or  even  a  fiber  with  a  strong  circular  birefringence  maintains  CV  vortex  at  large 
distance  [7]. 

Therefore  we  have  employed  either  fibers  with  a  slight  birefringence  or  circular  core- 
birefringent  ones  as  a  sensitive  element  of  a  fiber-optical  sensor  which  the  guided  vortices  can  be 
translated  through.  The  typical  scheme  of  a  OVFS  is  reproduced  in  Fig.l  [7]. 


Ill  is  the  vortex-forming  output  unit;  IV  is  the  detection  unit  matched  with  the  image  processing 
system.  1  is  the  He-Ne  laser;  2,17  are  the  beam  splitters;  3,6  are  the  mirrors,  4  is  the  A/4  plate;  5  is 
the  optical  wedge,  7,12  are  the  microlenses,  8,16  are  the  polarizers,  10  is  the  optical  fiber,  13,14- 
are  the  co-centered  microlenses,  15  is  the  collimator,  18  are  the  CCD  matrix  and  PC. 
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There  are  four  main  units  in  this  scheme  that  will  be  analyzed  in  the  following 
subsections:  the  fiber  excitation  unit,  the  vortex-guided  unit,  the  vortex-forming  output  unit  and 
the  detection  unit  matched  with  the  image  processing  system. 

The  fiber  excitation  unit  serves  for  the  excitation  of  the  CV  mode  in  a  low-mode  fiber. 

The  vortex-guided  unit  is  an  optical  fiber  positioned  into  a  special  arrangement 
reinforcing  influence  of  a  physical  parameter  under  test. 

The  output  unit  consists  from  two  microlenses  (and  diaphragm  between  them,  if 
necessary). 

The  detection  and  image-processing  unit  employs  the  CCD-matrix  matched  with 
computer  for  recording  and  processing  of  experimental  data  in  a  real  time. 

The  consequent  subsections  are  devoted  to  the  detailed  analysis  of  the  each  optical  unit. 
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4.2.2  High-efficiency  generation  of  optical  vortices 


As  a  rule,  coupling  of  an  optical  fiber  and  an  OV  beam  in  free  space  is  connected  with 
energy  losses  up  to  50%  and  more.  Besides,  energy  efficiency  of  an  OV  beam  generation  by 
computer-generated  holograms,  widely  discussed  in  a  number  of  works  (see,  for  example,  [1]  and 
references  therein),  amounts  several  percents.  Consequently,  a  computer-generated 
hologram/optical  fiber  system  has  very  low  efficiency  for  a  design  of  the  OVFS. 

On  the  other  hand,  a  number  of  the  experimental  results  [2]  shows  that  a  rather  simple 
optical  device  such  as  an  optical  wedge  enables  us  to  obtain  OVs  with  the  high  efficiency.  Some 
theoretical  aspects  of  this  problem  have  been  recently  considered  in  the  work  [3]. 

The  aim  of  the  given  section  is  the  experimental  and  theoretical  study  of  the  diffraction  of 
a  Gaussian  beam  by  an  optical  wedge  as  a  high  efficiency  device  for  an  OV  beam  generation. 

Let  us  consider  an  optical  wedge  illuminated  by  a  Gaussian  beam  (see  Fig.l)  in  such  a 
way  that  one-half  of  the  beam  passes  through  the  wedge  and  the  other-half  propagates  in  free 
space. 


7 


Fig.l.  The  Gaussian  beam  diffraction  by  the  optical  wedge. 

If  peripheral  rays  of  a  laser  beam  have  the  phase  difference  2k  ,  then  a  circuit  with  the 
radius  p  around  optical  axis  gives  2%  .  It  is  the  relevant  requirement  of  an  optical  vortex 

existence.  However,  the  phase  shift  on  the  computer-generated  hologram  is  not  dependent  on  the 
circuit  radius  while  this  requirement  does  not  fulfilled  in  our  case.  The  OV  formation  by  the 
optical  wedge  is  therefore  more  complicated  process. 

The  analysis  of  the  OV  generation  was  performed  by  means  of  the  Kirchhoff  integral  [4] 
for  a  far  field.  The  wave  function  has  a  form: 
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where  a  is  the  wedge  angle,  03  —  +  Zq  )>  Z0  —  — ~  is  the  Rayleigh  length, 

2  z 

erfc(z)  =  1-  erf(z),  erf{z)=  —r=  [ exp(-  / 2 )  dt  is  the  probability  integral. 

V7t  o 

The  key  role  in  the  OV  beam  formation  plays  the  relation  between  the  beam  waist  p  and  the 

wedge  angle  a.  From  eq.(l)  we  obtain  the  following  expression  for  a  very  wide  beam  (nearly 
plane  wave)  incident  on  the  wedge: 
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It  is  interesting  to  give  a  comparison  of  the  theoretical  results  obtained  by  means  of  eq. 
(2)  and  the  experiment  (see  below). 
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Fig.2.  The  generation  of  OVs  street  by  the  optical  wedge:  (a)  theory;  (b)  the  experimental  results 
obtained  by  a  microscope  cover  glass;  (c)  the  theoretical  phase  distribution  in  the  OVs  street.  The 
distance  from  the  glass  wedge  is  z  -  1000  X,  the  beam  waist  is  p  =  500  X ,  and  the  wedge  angle  is 
a  =  0.08  rad  . 


The  distance  Ay  between  neighboring  vortices  in  the  street  is  defined  by  Eq.  (2): 

X 

Ay  =  7 - itt - ,  (3) 

(n-l)tana 

where  n  is  the  refractive  index  of  the  wedge. 

A  single  OV  will  appearin  in  the  optimal  form,  if  the  distance  between  OVs  is  Ay  =  2p  . 


From  eq.  (3)  one  has 

X 

2(n-l)tana 


(4) 


The  good  agreement  of  analytical  (Fig.2a,c)  and  experimental  (Fig.2b)  results  allow  us  to 
rely  on  theoretical  predictions  in  designing  of  a  high  efficiency  device  of  the  OV  generation. 
Consider  now  the  properties  of  an  optimal-formed  OV. 
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Fig.3.  The  intensity  in  gray  scale  (a),  the  intensity  level  lines  (b),  phase  surface  (c),  and  equal 
phase  lines  (d)  for  the  optimal  OV,  calculated  by  eqs.(l)  and  (4)  under  z  —  1000A,,  a  =  0.08 rad 

The  experimental  investigation  of  the  OV  generation  by  an  optical  wedge  was 
implemented  in  experimental  set-up  shown  in  Fig.  4. 

The  first  set-up  (Fig.  4a)  enables  to  observe  the  generation  process  of  the  OV  street  and 
the  optical  wedge  features  by  means  of  the  interference  pattern  in  the  reflected  light  [5].  The 
second  set-up  (Fig.4b)  enables  to  study  the  properties  of  a  single  OV.  The  cover  microscope  glass 
plates  were  used  as  samples  of  optical  wedges.  The  interference  pattern  shown  in  Fig.5  points  out 
the  modulated  wedge  thickness. 


(a) 


(b) 


Fig.4.  The  experimental  set-up  for  the  investigation  of  a  Gaussian  beam  diffraction  by  an  optical 
wedge,  (a)  Generation  of  the  OV  street:  1  is  the  He-Ne  laser,  2  is  the  telescopic  objective,  3  is  the 
optical  wedge,  4,5  are  the  CCD  matrices,  (b)  Generation  of  a  single  OV :  1  is  the  He-Ne  laser,  2,8 
are  the  beamsplitters,  3,7  are  the  telescopic  objectives,  5,6  are  the  mirrors,  9  is  the  diaphragm,  10 
is  the  CCD-matrix. 
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(a)  (b) 


Fig.5.  The  interference  pattern  obtained  in  the  beam  reflected  by  the  cover  glass  (a)  and  the 
corresponding  OV  street  in  the  transmitted  beam  (b). 

Because  the  period  of  the  cover  glass  modulation  is  rather  large,  an  area  with  a  nearly 
straight  wedge  can  be  chosen. 

The  typical  experimental  picture  demonstrated  in  Fig.6  (a)  gives  a  good  agreement  with 
the  theoretical  prediction  of  a  single  OV  generation  depicted  by  Fig.6  (c). 


Fig.6.  The  intensity  distribution  and  level  lines  of  the  optimal  OV  generated  by  the  optical  wedge 
with  angle  a  =  8.8  X 10  ~fdd  :  (a,b)  the  experiment,  (c,d)  the  theory. 

The  considered  above  study  should  be  treated  as  a  preliminary  approach  to  the  problem  of 
a  high-efficiency  generation  of  OV  for  OVFS.  The  analytical  and  numerical  analysis,  confirmed 
by  the  experimental  results,  shows  the  simple  way  for  obtaining  of  an  OV  with  a  high  efficiency. 
The  shape  of  the  single  vortex  is  defined  by  both  the  wedge  and  laser  beam  parameters. 

References 

1.  M.S.  Soskin,  M.V.  Vasnetsov.  //  Horizons  in  World  Physics,  edited  by  M.  Vasnetsov  and 
K.  Staliunas,  1999.  V.  228.  P.  1-35. 

2.  N.  Alexeyev,  O.  S.  Borodavka,  V.  G.  Shvedov.  Birth  and  evolution  of  the  optical  dipoles  on 
the  dielectric  transparent  wedge,  Proceedings  of  SPIE.  1999,  3904.  P.68-73. 


11 


3 .  V.  Gorshkov  at  all.//  SPIE  proc.  2001 .  (to  be  published). 

4.  M.Bom,  E.  Wolf.  Principles  of  optics.  Oxford,  Pergamon  Press,  1968. 

5.  Izdebskaya  Ya.,  Shvedov  V.,  Kurabtzev  D.,  Alexeyev  A.,  Volyar  A.  “  The  optical 
vortex  generation  by  optical  wedge”  //  SPIE  proc.  (to  be  published). 

4.2.3  Design  and  assessments  of  the  two-arms  fiber-optical 
interferometer  on  the  base  of  an  optical  vortex 

The  main  advantages  and  imperfections  of  single-  and  multimode  sensors  having  been  in 
detail  discussed  in  the  Section  4.1.  The  aim  of  the  given  section  is  to  develop  the  design  and  to 
carry  out  the  assessments  of  the  OVFS  as  whole.  We  shall  consider  here  two  possible 
constructions  of  the  two-arms  interferometers  differed  by  transmission  ways  of  an  optical  vortex 
through  either  a  signal  or  a  reference  arm  [1],  In  the  former  case  the  OV  propagates  along  a  fiber 
and  it  is  a  guided  fiber  mode,  and  we  deal  with  a  free-space  OV  in  the  latter  case. 

The  greatest  advantage  of  the  OV  sensor  is  that  the  development  of  a  new  interferometric 
technique  of  the  image  processing  enables  us  to  extend  essentially  the  range  of  the  device 
linearity  and  to  eliminate  partially  noise  caused  by  intensity  instabilities  of  a  light  source. 

The  sketch  of  the  experimental  set-up  is  represented  in  Fig.  1 . 


Fig.l.  Sket(Tf^ 
the  OV  gen 
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laser,  S  Pr  is  the  beamsplitter,  T  is 
tie  mirror,  Sc  is  the  screen,  L  is  the  lens,  ^4 
MFisJhe  low-mode  fiber. 


The  linearly  polarized  light  from  the  single-frequency  laser  (k  =0.6328  pm.  power  is 
about  2  mW)  in  a  form  of  a  fundamental  Gaussian  mode  passed  through  the  beamspliter  forming 
the  two  arms  of  the  interferometer.  In  the  reference  arm,  the  beam  was  attenuated  by  the  polarizer 

P  and  then  transformed  into  the  circularly  polarized  one  after  the  %  -  plate.  The  reference  arm 

was  put  inside  a  thermostat  to  protect  the  arm  from  possible  ambient  temperature  variations.  In 
the  signal  arm,  three  main  optical  units  were  positioned.  In  the  vortex-generating  unit,  the 
Gaussian  beam  diffracted  by  the  optical  wedge  was  converted  into  the  OV.  The  scattered  light 
was  cut  off  by  the  optical  filter  consisting  from  two  lenses  and  a  diaphragm  between  them. 

The  cleaned  OV  was  launched  into  the  fiber  by  the  20x  microlens.  The  conventional 
low-birefringent  optical  fiber  consisting  from  the  G  e02 1  Si02  core  and  the  B203ISi02 
cladding  with  circular  cross-section  was  used  as  a  guided  vortex  system.  The  birefringence  was 
less  than  icr6  ,  the  core  radius  was  p  =  4  |Lim,  the  fiber  parameter  was  V  =  5.1.  This  fiber  could 
maintain  the  following  modes  (in  the  linearly  polarized  basis):  fundamental  H  ,  H  E°±  and 
HEi2,HE°2  modes,  azimuth  symmetric  TE01,TM01  modes  and  hybrid 
E  H  *i,E  H  ,H  E  zlfH  E  ^  modes.  They  are  associated  with  the  following  types  of  vortices: 

HE  n  ,H  E^,CV  1\,  IV n,C  V2l .  The  higher  modes,  such  as  the  CV2 L  vortex  or 
E  H 11  ,E  H  ,H  E  ,H E  31  modes,  have  the  excitation  effency  less  than  10  %.  The  fiber  ends 
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were  fixed  on  a  3 -coordinate  tables  provided  by  microscrews  with  transverse  and  longitudinal 
displacements  about  5  pm.  The  input  and  output  microobjectives  were  mounted  on  the  3- 
coordinate  tables  provided  by  microscrews  and  a  hydraulic  mechanism  with  accuracy  up  to  0.5 
pm.  The  OV  generation  unit  has  also  the  lateral  image  microscope  to  observe  the  vortex  spot 
position  at  the  fiber  input  face.  The  fiber  was  coiled  without  tension  around  a  copper  cylinder 
with  32  mm  outer  diameter  and  30  mm  inner  diameter  was  placed  into  a  bath  with  oil.  Then,  this 
construction  was  put  in  a  thermostat,  where  the  temperature  was  measured  by  a  thermocouple 
with  the  precision  0.05°  C  at  the  range  from  20°  to  50°  C. 

Two  centered  20x  microobjectives,  mounted  on  optical  tables,  were  placed  at  the  fiber 
output.  The  optical  tables  were  supplied  with  microscrews  and  a  hydraulic  arrangement  to  provide 
an  accuracy  of  the  displacement  up  to  0.5  pm. 

The  signal  and  reference  beams  were  combined  by  the  beamsplitter  and  directed  to  the 
window  of  a  CCD-matrix  producing  the  image  on  the  TV  monitor.  The  interference  pattern  looks 
like  a  spiral  fringe  [1],  The  CCD-matrix  having  512x582  sensitive  elements  and  was  matched 
with  computer  Pentium-3  (800)  to  process  the  spiral  interference  image  in  real  time.  We  detected 
the  rotation  of  the  interference  spiral  and  calibrated  the  rotation  angle  in  temperature  degrees. 

The  experiment  runs  as  follows.  The  OV  generation  unit  forms  a  CV  mode  for  launching 
it  into  the  fiber  (see  Section  4.2.3).  The  free  OV  and  the  guided  CV  mode  will  be  completely 
coupled  if  only  their  transversal  parameters  are  the  same.  In  reality  the  OV  beam  axis  is  slightly 
inclined  and  a  little  bit  displaced  relative  to  the  fiber  axis.  As  a  result,  the  additional  eigen-modes 
in  the  fiber  such  as  IV  mode  and  the  fundamental  H  En  mode  are  excited.  The  H  En  mode  field 

displaces  merely  the  CV  mode  center  and  causes  no  noise  for  the  sensor  operation.  The  IV  mode 
produces  fairly  strong  noise  effect  to  the  interference  image.  Even  a  slight  displacement  of  the 
OV  beam  spot  relative  to  the  fiber  axis  (less  than  0.5  pm)  stimulates  an  appearance  of  the  IV 
mode  (see  Section  4.2.4).  The  sensitivity  of  the  sensor  radically  reduces  as  the  measurement  error 
raises.  Besides,  even  if  the  free  and  guided  vortices  are  strongly  coupled  and  the  only  CV  mode  is 
excited  in  the  fiber  input,  any  slight  perturbations  of  the  fiber  itself  cause  the  IV  mode  birth  and, 
consequently,  growth  of  the  measurement  error.  For  the  pure  CV  mode  exitation  the  laser  spot  at 
the  fiber  face  must  be  positioned  with  accuracy  better  than  0.5  pm  relative  to  the  fiber  axis. 

If  the  only  CV  mode  appears  at  the  fiber  output,  or  the  noise  intensity  is  as  little  as  third 
of  the  CV  mode  intensity,  the  filtering  two-lens  system  is  removed  from  the  set-up.  For  large 
noise  level  the  measurements  are  run  as  it  is  described  in  Section  4.2.4. 

After  adjusting  the  set-up  we  carry  out  cleaning  of  an  interference  spiral  image  from  the 
noise.  The  contrast  of  the  spiral  image  can  be  improved  by  equalizing  of  the  intensities  of  the 
reference  and  signal  beams.  Besides,  the  interference  image  must  have  no  large-scale  patterns.  To 
eliminate  them,  optical  filter  is  placed  in  the  reference  arm.  In  the  signal  arm  this  function  carries 
out  the  two-lenses/diaphragm  system. 

The  image  of  the  interference  spiral  is  directed  to  the  computer  image  processing.  The 
experimental  results  (temperature)  are  displayed  on  the  TV  monitor  simultaneously  with  the  spiral 
image  in  a  real-time  scale  with  the  time  resolution  about  0.4  c  (see  Section  4.2.5). 

Figure.2  represents  the  typical  experimental  results. 
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Fig.2.  The  experimental  dependence  of  the  spiral  angular  position  (p  (in  rad)  vs  the  temperature 
t°  C  for  the  fiber-optical  interferometer  with  the  CV  in  the  signal  arm. 


As  it  can  be  seen  from  the  curve  in  Fig.  2,  the  temperature  dependence  of  the  angular 
rotation  has  a  linear  character  at  a  temperature  range  from  20  to  50  0  C  .  In  contrast  to 
conventional  types  of  fiber  optical  sensors,  our  vortex  sensor  has  a  linear  range  of  the  angular 
dependency  within  0-r  2n  interval  of  phase  shift  variations  between  a  referent  and  signal  beam. 


Deviations  of  the  linearity  are  within  to  0.01  rad.  The  temperature  sensitivity  is  about  140  - . 

Km 

The  experimental  reproducibility  for  different  2.5-meter  pieces  of  the  same  fiber  is  up  to  2  %. 

We  have  also  considered  an  alternate  type  of  a  OVFS  [2].  A  single-mode  fiber  was  used 
in  a  signal  arm  and  a  free  OV  beam  served  as  a  reference.  The  single-mode  low-birefringence 
fiber  had  the  following  parameters:  the  core  radius  p  =  3.4  pm,  the  waveguide  parameter  V=  2.1 

at  the  wavelength  0.63  pm,  the  linear  birefringence  was  less  than  10  6 .  So,  the  fiber  can  maintain 


the  HE ^  mode  only.  For  example,  exciting  the  fiber  by  the  right-hand  polarized  with  an 
ellipticity  degree  Q  =  0.99  Gaussian  beam,  we  obtain  Q  =  0.73  after  passing  through  the  10  m- 
length  fiber.  Consequently,  the  beating  length  in  the  given  low-birefringent  waveguide  is  less  than 
20  m. 
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Fig. 3.  The  experimental  curve  as  a  function  of  the  interference  spiral  rotation  cp  (in  rad)  vs  the 

temperature  t 0  C  for  the  fiber-optical  interferometer  with  a  vortex  in  a  reference  arm;  fiber  length 
is  about  27  cm. 

The  free  OV  propagated  along  the  signal  arm.  It  should  be  noted  that  a  coherence  length 
of  our  single-frequency  laser  source  was  a  fairly  large.  The  contrast  of  the  interference  image  in 
the  fiber-optical  interferometer  with  the  arm  difference  about  30  m  was  not  less  than  0.9. 

Similar  to  the  former  case  we  measured  the  dependence  of  the  angular  rotation  of  spiral 
image  vs  temperature  variations.  Figure.3  illustrates  a  nearly  straight  line  as  a  function  of  the 
angular  position  (p  (in  rad)  and  the  temperature  t°  C.  Though  the  latter  design  is  a  steady  sensor 

rad 

relative  to  random  mechanical  perturbations  we  can  see  that  the  sensibility  of  it  is  about  110  - 

Km 

while  the  linearity  and  the  accuracy  is  similar  to  the  former  sensor. 

The  grater  sensibility  of  the  interferometer  with  a  vortex  in  a  signal  arm  is  a  result  of  that 
the  guided  CV  vortex  propagates  along  a  fiber  with  a  less  phase  velocity  than  a  H  Enthat  is  an 
optical  path  for  the  same  fiber  length  is  larger. 

Thus,  we  have  experimentally  realized  two  types  of  a  two-arm  OVFS  with  the  OV  in  the 
signal  or  reference  arm  and  carried  out  the  assessments  of  their  operation. 
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4.2.4  Fiber-optical  interferometer  with  the  OV  recovery  system 

Notwithstanding  the  fact  that  any  slight  ambient  perturbations  affected  the  OV  on  the 
base  a  conventional  low-mode  fiber  cause  a  strong  mode  noise  of  the  interferometer,  its  high 
sensitivity  and  a  constructive  simplicity  permit  to  design  special  devices  for  filtering  of  outside 
noise.  At  the  same  time,  in  order  to  avoid  of  unnecessary  energy  losses,  a  noise  field  must  be 
conversed  into  a  useful  signal.  Evidently,  there  exist  no  physical  mechanisms  to  ensure  that 
process  over  all  optical  fiber  length.  Nevertheless,  it  had  been  known  some  nonlinear  ways  for 
recovering  of  a  useful  optical  information  without  energy  losses  at  the  fiber  input  [1],  In  our  case, 
it  is  necessary  to  recover  an  optical  information  carrying  with  optical  vortices.  The  fact  is  that  an 
initial  guided  OV  split  into  a  number  of  other  optical  vortices  at  the  expense  of  external 
perturbations  but  a  total  topological  charge  is  preserved.  It  turns  out  to  be  possible  an  inverse 
optical  process  of  annihilating  of  «parasitic»  optical  vortices  and  recovering  of  the  initial  one  so 
as  a  total  topological  charge  is  preserved  as  well  [2,3]. 

The  aim  of  the  given  section  is  to  study  experimentally  and  theoretically  death,  birth  and 
recovery  processes  of  optical  vortices  at  an  optical  fiber  input. 

Let  us  consider  the  optical  system  consisting  of  an  optical  fiber  and  two  short-focusing 
lenses  (Fig.l).  The  first  lens  forms  the  image  of  the  fiber  output  at  the  inner  observation  plane  Z, 
the  second  lens  passes  this  image  at  the  outer  observation  plane  P. 


z 
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Fig.  1. The  sketch  representation  of  the  experimental  set-up:  lis  the  optical  fiber,  2  are  the  short- 
focusing  lenses 

Let  us  suggest  that  the  CV  vortex,  the  fundamental  H  En,  H  Eu  and  LP22  modes  are 

realized  in  the  optical  fiber.  It  is  a  typical  case  of  the  CV  excitation  in  a  conventional  low-mode 
fiber  [4,5].  In  the  scalar  nonparaxial  approximation  the  wave  function  of  the  light  near  the  focal 
plane  Z  has  a  form  [2,3]: 


vF(d)  =  ^“4  G  •  'F1(f>) + F  •  VP(J0)  +D •  44 


(1) 


where  44  = 


1  ~\kz{)  \i<MR)-}-kz( 


r2  ^ 

2-3  — 
R2 


j2(kR )  is  the  combined  Legandre-Bessel 


beam,  44  = — coscp  /',  (kR)  is  the  Legendre-Bessel  beam  with  l  -  m  -  1,  4°  *  -  jo(kR)  is 
R 


the  fundamental  Bessel  beam,  — exp(/(p)  /,  (kR)  is  the  nonparaxial  OV,  jm{kR )  is  the 

R 

spherical  Bessel  function  of  the  first  kind,  m-th  order,  R  =  -Jr2  +  (z  +  )2  is  the  complex 


radius,  G,F,D  are  the  constant  coefficients. 

Two  possible  wave  processes  take  place  in  the  system  in  accordance  with  the  value  of  the 
coefficient  D  in  eq.  (1).  If  D=0,  the  fiber  is  excited  by  a  topologically  neutral  beam  and 
symmetric  topological  dipoles  or  quadruples  [6]  appear  at  the  fiber  output.  (The  symmetric 
dipoles  or  quadruples  represent  a  system  with  the  same  «weights»  of  positive  and  negative 
topological  charges.)  If  D  ^  0 ,  there  are  a  topological  charged  beam  at  the  fiber  input. 

The  Fig.  1  shows  the  case  of  the  topological  neutral  fiber  excitation.  The  space  trajectories 

were  built  from  the  condition:  R^'P(x.y.z)]  =  0 ,  I  m[vP(x,  y,  z)]  =  0  and  characterize  the 
behavior  of  the  symmetric  topological  quadruple. 


(b) 


Fig.2.  The  space  trajectories  of  the  symmetric  topological  quadruple  (D=0):  (a)  F=l,  G=3;  (b)  F=- 

1,G=2. 
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The  nonparaxial  topological  quadruple  has  two  types  of  trajectories.  They  differ  by  forms 
of  forbidden  or  allowed  zones.  The  trajectory  in  Fig. 2a  corresponds  to  the  quadruple  within 
allowed  zones.  At  the  front  wall  of  the  zone,  four  coupled  OVs  of  the  quadruple  are  bom  passing 
through  the  focal  plane.  They  annihilate  at  the  back  wall.  There  are  no  OVs  outside  of  the  zone. 

The  opposite  case  is  shown  by  the  space  trajectory  in  the  Fig.2b.  There  is  the  forbidden 
zone,  at  the  edges  of  which  the  OVs  are  annihilated  or  bom.  There  are  no  OVs  inside  of  the  zone. 


Fig.3.  The  asymmetric  topological  quadmple:  F=0.5,  G=-0.3,  D=l. 


Fig.4.The  space  trajectories  of  the  asymmetric  dipole:  (a)  the  hyperbolic-like  trajectory  for  F=l, 
G=-l,  D=l;  (b)  the  parabolic-like  trajectory  for  F=-0.5,  G=-0.1,  D=l;  (c)  the  elliptic-like 
trajectory  for  F=-l,  G=0.1,  D=l. 

The  practical  interest  presents  the  case  of  a  topological  charged  excitation  beam  (D  ^  0 ). 
The  space  curves  of  the  Fig.3  represent  the  OV  behavior  of  the  beam  carrying  both  coupled  and 
free  vortices.  Similar  to  the  former  case  there  are  both  forbidden  and  allowed  zones  but  they  exist 
simultaneously.  Besides,  the  quadruple  becomes  the  asymmetric  one.  At  the  front  edge  of  the 
forbidden-allowed  zone,  the  four  OVs  are  at  first  annihilated  and  then  bom. 

However  «the  weights»  of  these  OVs  are  unequal.  The  OV  with  the  topological  charge 
identical  to  the  initial  OV  is  transmitted  near  the  optical  axis,  whence  the  rest  OVs  are  forced  out 
far  from  the  optical  axis  and  have  a  negligibly  small  «weight»  in  the  vicinity  of  the  focal  plane. 
At  the  back  edge  of  the  forbidden-allowed  zone  the  OVs  are  at  first  annihilated  and  then  bom.  As 
the  singular  beam  is  transmited  far  from  the  focal  plane,  the  «weights»  of  OVs  are  equalized. 
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Fig.4  represents  other  calculated  practically  interesting  situations.  It  is  important  to  keep 
in  mind  that  the  OV  in  the  asymmetric  dipole  and  quadruple  is  selected  near  the  focal 
plane.Therefore,  the  rest  coupled  OVs  are  either  annihilated  or  their  «weight»  is  essentially 
reduced,  and  they  are  forced  out  far  from  the  optical  axis.  In  any  single  case,  the  «parasitic»  OVs 
are  suppressed,  and  information  is  carried  by  the  useful  OV  without  energy  losses.  The  numerical 
estimations  show  that  the“distance  of  life”  of  the  free  OV  is  about  some  wavelengths  at  the 
vicinity  of  the  beam  waist. 

As  a  base  of  the  experimental  set-up,  the  device,  schematically  represented  in  the  Fig.l 
from  .4.2.2,  was  used.  The  main  unit  of  the  device  was  the  centered  system  of  two  40x 
microobjectives  fixed  on  the  three-coordinate  tables.  The  microsrew  system  and  the  hydraulic 
arrangement  permit  us  to  displace  the  microobjective  along  the  optical  axis  and  in  the  transverse 
directions  to  within  0.5  pm.  As  before,  the  conventional  fiber  was  coiled  on  a  cylinder  and 
excited  by  a  circularly  polarized  light  from  the  He-Ne  laser.  Beforehand,  the  polarized  Gaussian 
beam  passes  through  the  optical  wedge  and  converts  into  an  OV.  Besides,  the  signs  of  the 

topological  charge  and  the  circular  polarization  form  either  |  +  l,+l)  or  group.  Such 

adjustment  of  the  initial  beam  ensures  the  primary  excitation  of  the  CV  in  the  fiber.  In  order  to 
insure  the  topologically  neutral  fiber  excitation,  it  is  sufficiently  to  remove  the  optical-wedge 
transparence  from  the  beam. 

Concerning  the  output  radiation,  the  analyzer  may  be  inserted  in  the  beam  in  order  to 
make  better  an  image  contrast.  The  fundamental  Gaussian  beam  in  the  reference  arm  has  a 
circular  polarization.  The  CCD-matrix  detector  has  512x582  sensible  elements  recording  both  the 
signal  beam  and  the  interference  pattern.  The  experimental  results  are  operated  by  the  computer- 
Pentium  3(800)  in  real-time  scale. 


Z,mcm 
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Fig.5.  Space  trajectories  of  the  coupled  and  free  OVs  in  the  output  beam  near  the  focal  plane: 

(a)  the  topologically  neutral  initials  beam;  (b)  the  topologically  charged  initials  beam. 

The  radiation  after  the  optical  fiber  turns  out  to  manifest  itself  in  special  way  unexpected 
before.  It  had  been  predicted  that  the  beam  after  the  fiber  matches  completely  from  the  Laguerre- 
Gaussian  modes  of  free  space  and  propagates  to  the  far  field  without  any  structural  changes.  Our 
experiments  show  that  such  is  not  the  case  [7], 

We  have  investigated  the  evolution  the  beam  intensity  distribution  with  the  distance 
between  the  observation  plane  Z  (see  Fig.  1)  and  the  end.of  the  fiber  for  two  cases,  when  the 
incident  beam  is  topologically  neutral  or  charged.  The  former  the  latter  casesaredepicted  in  the 
Fig.5a,b.  The  upper  pictures  are  the  interference  patterns  with  typical  fork’s  positions.  It  is  those 
positions  that  represent  the  bellow  curves  as  the  functions  of  the  distance  R  from  the  optical  axis 
to  the  fork  and  the  distance  z  between  the  fiber  end  and  the  observation  plane. 
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Fig.6.  The  interference  patterns  of  the  symmetric  (a)  and  asymmetric  topological  quadruples  and 
the  recovered  OV  (c).  For  comparison  the  interference  spiral  of  the  “ideal”  computer-generated 
OV  is  depicted  (d). 

The  space  curve  in  the  Eg.5a  characterizes  the  positions  of  the  coupled  OVs  in  the 
topologically  neutral  beam,  The  fiber  excitation  is  chosen  by  such  a  way  that  the  coupled  OVs 
form  the  elliptic-like  trajectory  similar  to  that  depicted  in  Fig.4c.  The  form  of  the  trajectory  points 
out  existence  of  the  forbidden  zone  with  extension  about  50  pm.  At  its  edges  OVs  are  born  and 
annihilated. 

A  completely  different  type  of  the  space  trajectory  occurs  in  the  topologically  charged 
beam  (see  Fig.5b).  Indeed,  there  are  two  closed  pieces  of  the  curve  positioned  before  and  after  the 
focal  plane  and  one  straight-line  crossing  the  beam  waist.  The  closed  curves  are  produced  by  the 
coupled  OVs  while  the  straight  line  is  traced  by  the  free  OV.  The  free  OV  is  located  inside  the 
forbidden  zone  with  “life  length”  about  3  pm.  It  enables  to  ensue  suppressing  of  the  mode  noise 
in  a  OVFS  without  energy  losses. 

Although  the  intensity  distribution  of  the  recovered  field  has  not  an  ideal  structure  of  the 
OV,  the  interference  patterns  are  near  to  spiral  of  the  ideal  OV.  Fig.6  illustrates  a  typical 
interference  image  of  the  recovered  field. 

Thus,  we  have  shown  that  the  centered  microobjective  system  placed  at  the  fiber  output  is 
able  to  recover  the  image  of  a  single  OV  without  energy  losses.  Similarity  of  the  recovered 
interference  spiral  and  the  ideal  one  permits  to  design  operating  samples  of  the  OVFS  on  the  base 
of  conventional  commercial  optical  fibers. 
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4.2.5  Computer  processing  of  the  interference  spiral 

A  computer  processing  of  an  operative  information  from  fiber-optical  sensors,  as  a  rule, 
is  succeeded  by  operation  with  numerical  data,  for  example,  a  total  intensity  after  a  sensor  [1,2]. 
At  the  same  time,  it  is  well  known  that  intensity  fluctuations  of  a  laser  source  and  other  similar 
reasons  cause  the  greatest  errors  in  operation  of  a  device  reducing  both  accuracy  and  sensitivity. 
On  the  other  hand,  a  destructive  effect  of  the  initial  intensity  variations  can  be  essentially  reduced 
or  entirely  eliminated  by  means  of  processing  of  an  interference  image  as  a  whole,  by  a  computer 
image  processing. 


20 


The  aim  of  the  given  section  is  to  develop  and  to  test  the  computer  image  processing 
technique  of  the  angular  location  of  the  spiral  interference  pattern  as  a  whole  in  a  vortex  fiber- 
optical  sensor. 

The  phase  shift  between  singular  and  reference  beams  in  the  ideal  case  is  a  linear  function 
of  the  angular  position  of  the  interference  spiral.  Consequently,  optical  fiber  perturbations  caused 
by  external  parametric  fields  (a  temperature,  a  pressure,  acceleration,  etc.)  can  be  defined  by 
means  of  the  direct  measurement  of  the  spiral  angular  position.  Moreover,  in  contrast  to  other 
interference  ways  this  one  gives  more  high  operating  accuracy  and  reproducibility  of 
experimental  results. 

Consider  the  intensity  distribution  occurring  as  a  result  of  the  interference 
between  the  optical  vortex  and  Gaussian  beam  as  [3] : 

I  =  1 1  +  1 2  +  cos($!  —  ^2)  ’  (1) 

where  I,,  I2  <f>2  are  the  Gaussian  beam  and  optical  vortex  intensities  and  phases,  respectively. 
The  interference  image  looks  as  bright  and  dark  spirals.  Position  of  the  bright  curve  is  defined  by 
the  phase  difference: 


<P1 -<1>2  =  2mn,  (m  =  0,1,2...) 

(for  the  dark  curve,  it  is  O,  —  ,  =  (2m  +  1  )7t ).  The  Gaussian  beam  phase  is  [5] 
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length,  is  the  beam  waist,  is  the  initial  phase. 
On  the  other  hand,  the  OV  phase  is  [4] 
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is  the  wavefront  radius  curvature,  z0g  —  is  the  Rayleigh 
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where  Rv\z)  —  z  1 H — j-  is  the  OV  curvature,  z0v  =  — —  is  the  Rayleigh  length,  pv  and  0,,,  are 

V  Z  J  2 

the  OV  waist  radius  and  initial  phase,  r,  9  are  the  polar  coordinates  in  the  observation  plane. 
From  equations  (2)  -  (4)  we  have 
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Since  z-const  at  the  observation  plane,  the  equation  (5)  can  be  rewritten  as: 
<3>j  — <E>2  =-ar2  +b'  +  /cp  , 
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The  equation  of  equiphase  lines  on  the  spiral  image  can  be  written  as  lq>  =  ar 2  +  b0 

(b0  —  CP .  —  <J>1  +b').  The  most  simple  spiral  form  being  comfortable  for  the  image  computer 
processing  (/=/)  expressed  as: 


9  =  ar2  +  b0 . 


(V) 


We  have  considered  various  techniques  of  the  computer  processing  of  spiral  images. 
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The  first  technique  defines  the  phase  difference  as  the  position  (r,  cp )  of  the  intensity 
minimun  along  circumference  with  the  given  radius  (see  fig.l). 


The  intensity 
distribution 
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the  intensity 
minimum 


r  i 
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Fig.l.  The  first  technique  of  the  spiral  image  processing. 


The  advantage  of  this  technique  as  against  the  other  ways  is  the  very  high  rate  of  the 
experimental  data  processing,  while  its  imperfection  is  the  high  sensitivity  to  a  local  noise. 

The  second  technique  defines  the  angular  position  of  the  line  with  a  major  minimum 
intensity  obtained  as  a  Moire’s  pattern. 

The  Moire’s  patterns  are  obtained  by  subtraction  of  two  interference  images:  initial  and 
instantaneous.  The  inversion  process  is  operated  relative  to  the  X-axis  of  the  mirror-reflection  of 
the  image. 

The  obtained  image  has  the  strongly  delineated  minimum  in  a  form  of  the  black  band  (see 
Fig.2).  The  theoretical  value  of  this  minimum  can  be  defined  as  follows.  Due  to  a  mirror- 
reflection,  the  angle  (p  transforms  into  -cp  in  eq.  (6).  It  corresponds  to  the  sign  change  of  the  OV 
topological  charge  and  the  sign  of  the  curvature  of  the  spiral  changes  too.  Mathematically  this 
process  is  written  as: 

/  =  7j  +  1 2  T  2^V7 cos(ar2  —  /cp  +  ([) )  —  {/j  +  /2  +  2^JT^I~2  cos (c?r2  +  /cp J  (8) 


(/=1)  or 

I  —  2-n//1/2  {cos(ar2  -  /cp  +  c|) )  —  cos(ar2  +  /cp)}=  4^7, /2  sin|^ar  jsin{/cp  —  c|>  /  2} 

d) 

As  it  can  be  seen  from  the  last  expression,  an  axis  of  the  black  band  turns  into  the  —  angle 

2 

relative  to  X-axis. 


22 


The  tnU]  intensity 
aiong  □  radius 

1 

7 

^ - J 

The 

miminiTin 

azimuth 

coordinate 


Fig.2.  The  second  technique  of  the  image  processing. 

The  advantage  of  the  given  technique  as  against  the  first  one  is  the  noise  diminution  at  the 
expense  of  the  intensity  averaging-out  over  the  radius.  Its  imperfection  is  the  large  sensibility  to 
any  shifts  of  the  spiral  center.  Deviations  of  the  center’s  coordinates  cause  appearance  the  “fork¬ 
like”  intensity  distribution  instead  of  the  clearly  delineated  straight  band  (see  Fig.3). 


Fig.3.  The  fork-like  image 


The  third  technique  defines  the  phase  position  cp0  along  the  darkest  line  in  the 
interference  image  (see.  Fig.4). 

As  the  matter  of  fact,  the  interference  image  is  constellation  of  the  dark  and  the  bright 
spirals.  The  brightness  of  these  spirals  is  defined  by  the  phase  difference  between  beams  on  the 
given  curve.  Knowing  the  spiral  curvature  parameter  a,  we  can  always  define  the  integral 
intensity  on  the  given  spiral  and,  consequently,  obtain  the  required  phase  shift.  This  technique 
underlains  the  experimental  measurement  of  the  real  phase  differences  induced  by  the  external 
parametric  fields. 

The  main  advantage  of  the  given  technique  is  the  very  weak  sensitivity  to  ambient  noise 
due  to  the  intensity  integration  along  the  spiral  curve.  Besides,  it  has  the  lesser  sensitivity  to  the 
determination  of  the  interference  spiral  center  as  against  the  2-nd  way. 

In  order  that  these  techniques  can  be  realized  it  is  necessary  to  find  coordinates  of  the 
interference  spiral  center.  Find  “the  center  of  gravity”  of  the  dark  and  the  bright  spirals.  They  are 
a  little  bit  removed  relative  to  the  center  of  a  true  image.  But  “the  center  of  gravity”  of  the  dual 
spiral  is  just  the  sought  true  center  of  the  all  interference  pattern. 
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Fig.4.  The  third  technique  of  the  image  processing. 


For  that  reasons  the  next  mathematical  method  was  developed.  The  interference  image 
was  scanned  along  the  X-  or  Y-  axes.  On  the  each  line,  the  intensity  minimum  and  maximum 
were  defined.  Besides  the  xi  -  and  yi  -  coordinates  of  the  midst  of  the  most  wide  interval  between 
the  maximum  and  minimum  was  found.  Then  the  coordinates  the  spiral  center  were  expressed  as 

N  N 


X.  - 


F  - 


(9) 


N  N 

where  N  is  the  number  of  lines  with  the  most  wide  interval. 

The  x,  coordinate  of  the  widest  interval  is  a  little  removed  relative  to  the  true  center. 


However  there  is  any  other  mirror-symmetric  line  relative  to  the  true  center.  Moreover  its  x' 
coordinate  on  this  interval  is  removed  at  the  same  distance  but  with  opposite  sign.  The 

X  +  X 

superposition  — - gives  the  location  of  the  true  spiral  center  (see  Fig.5). 


xp — IDO-Vr) — 1 72 


Fig.5.  The  detecting  of  the  spiral  center  coordinates. 


To  realize  the  third  technique  [6]  it  is  necessary  to  find  the  a-  parameter  of  the  spiral 
curvature  from  the  eq.(7).  To  do  that,  it  can  be  used  the  method  similar  to  the  1-st  one.  Let  us 
change  the  radius  of  some  ring  encircling  the  spiral  center,  from  rmin  to  rmax .  we  find  the 
intensity  minimum  with  the  azimuth  angle  (p  on  the  each  circle.  The  Fig.  6  shows  that 
dependence  cp  (r2 ). 
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Fig.6.  The  detecting  of  the  a-  parameter  of  the  spiral  curvature. 


The  c/ -parameter  of  the  spiral  curvature  was  found  by  means  of  the  least-squares  method. 
The  a-  parameter  and  the  spiral  curvature  K  are  coupled  by  each  other  as 


r[4(cp— cp0)2+3l/4(cp-cp0) 

v[4(<P-<Po)  +  1l3 


The  experimental  data  as  functions  of  temperature  t"C  was  obtained  on  the  base  of 
OVFS.  The  CCD-camera  fixs  the  changes  of  the  interference  image.  The  temperature  varies  from 
18.4  ‘N  up  to  24.6  ‘N  during  10  min.  The  temperature  variations  from  one  frame  to  another  was 
AT~4  10"4 '  N.  The  dependence  of  the  spiral  phase  on  the  number  of  frame  (it  corresponds  to  the 
temperature  variations)  for  the  frame  succession  was  computed  in  concordance  with  the 
techniques  1-3.  The  Fig.7  represents  these  results  for  the  techniques  1-3.  The  obtained 
dependence  was  approximated  by  the  linear  function  y-ax+b.  Further,  the  standard  deviation  was 


defined  as 


Py,~y) 

|  n(n- 1) 


.  Values  of  the  parameters  a,  b,  A  are  represented  below  the 


curves. 

The  comparison  of  these  curves  shows  that  the  best  accuracy,  linearity  and  also  the  least  deviation 
corresponds  to  the  third  method. 

To  process  the  interference  image  in  a  real-time  scale,  we  have  developed  the  specific 
computer  program  adopted  to  DELPHI-5,  what  enables  us  to  realize  the  image  processing  of  one 
frame  during  0.4  c  in  according  to  the  third  technique  (see  Fig.8).  For  the  beginning,  the  program 
defines  the  coordinates  of  the  spiral  center  and  the  curvature  parameter  a  averaged  over  10  frames 
(see  Fig.8a).  Provided  their  mean  square  is  more  than  10%,  the  process  is  repeated  until  these 
values  will  be  10%  or  less  to  eliminate  ambient  noise  hindrances  exceeding  a  given  threshold. 
Further,  an  angular  position  of  the  spiral  is  viewed  on  the  TV  monitor  graduated  in  terms  of 
temperature  (see  Fig.8b). 
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Fig.7.  The  dependence  of  the  spiral  phase  on  the  number  of  frames 
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Fig.8.  View  of  the  operating  windows  of  the  computer  processing:  (a)  processing  of  the  spiral’s 
parameters,  (b)  processing  of  the  image  angular  position 

Thus,  we  have  developed  and  compared  three  techniques  of  the  computer  image 
processing  for  the  OVFS.  It  was  shown  that  the  third  technique,  i.e.  determination  of  the  angular 
position  along  the  dark  line  in  the  interference  image  has  been  employed  as  most  optimal.  The 
developed  technique  permits  to  process  a  spiral  interference  image  in  a  real-time  scale  during  0.4 
c  per  one  frame.  The  accuracy  of  the  temperature  measurement  is  up  to  0.012  rad  in  the  angular 
scale. 
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4.2.6  Design  and  assessment  of  one-arm  OVFS 


The  main  problems  with  two-arm  interferometers  such  as,  for  example,  a  single-mode 
Mach-Zehnder  interferometer  configuration,  are  attributed  to  the  noise,  especially  instability  and 
spurious  phase  shifts.  The  one-arm  interferometric  configuration  holds  undoubted  advantages  in 
this  point.  The  fact  is  that  the  interfering  modes  transmitting  through  the  same  fiber  possess 
exactly  the  same  ambient  perturbations.  The  only  distinction  appears  due  to  a  different  response 
of  the  mode  propagation  constants  to  temperature  deviations. 

The  aim  of  the  given  section  is  to  study  theoretically  and  experimentally  the  optical 
properties  of  the  one-arm  OVFS  on  the  base  of  the  guided  OV  and  to  carry  out  its  assessment. 

As  it  was  mentioned  above,  the  main  guided  modes  in  a  fiber  are  CV,  IV  and  H  E  wave 
fields  with  circular  polarization  [1,2].  In  contrast  to  the  CV  mode,  IV  one  manifests  itself  as  the 
wave  superposition  of  the  azimuthall-symmetric  TE01  and  TM  01  modes.  The  “weights”  of  these 
modes  are  defined  by  the  excitation  conditions  and  they  can  be  strongly  specified  at  the  fiber 
input.  The  IV  mode  being  broken  up  into  the  eigen-modes  may  change  both  the  topological 
charge  and  the  sign  of  its  circular  polarization. 

Consequently,  the  circular  polarization  analyzer  at  the  fiber  output  can  suppress  the  CV 

mode  and  H  E^  mode,  but  not  the  IV  mode. 

Let  us  suggest  that  the  H  E  ^  mode  and  CV,  IV  modes  are  simultaneousely  excited  at 
the  fiber  input.  The  wave  function  has  a  form  0  =>  A  ■  CVf  +  B IV f  +  C  HEf  or 
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where  Fl  0  (R)  is  the  function  of  the  field  profile,  R  =  — ,  (3  is  the  scalar  propagation  constant 

p 

g|^  _ gP 

of  the  CV  or  IV  modes,  8(3  (T  is  the  polarisation  correction  to  CV  mode,  A (3  =  — - — , 

5[3/u  TE  are  the  polarisation  corrections  for  TE  and  TM  modes,  (3/f/  is  the  propagation  constant 

of  HEU  mode,  A,  B  and  C  are  constant  coefficients. 

If  the  beam  after  the  fiber  and  and  circular  analyser  loses  the  left-hand  polarisation,  the 
zero-intensity  point  moves  according  eq.(l): 


{exp(/cp)+  Gexp(-/cp)cosA|3  exp(/5|3  z)}  Fl(R)+  Fexp(/8(30z)F0(i?)  =  0, 


(2) 


B 


C 


where  8|30  =(3ff£-(3  -5(3CT,5(3  =A(3-8|3cy,  G  =  —,  F  =  —. 

Let  us  consider  three  possible  cases. 

1)  The  CV  and  H  E  ^  modes  are  excited  in  the  fiber  only.  Then 


cp  =7t  +5(30z,  Fl(R)+FFn(R).  (3) 

It  means  that  the  OV  center  traces  a  circular  line.  The  angular  position  of  the  amplitude  zero 
depends  both  on  the  distance  z  and  the  difference  8(30  between  the  propagation  constants.  The 
temperature  sensitivity  of  the  angle  position  of  the  zero  point  is 

3- 


or 


Acp  3(5(3  0 ) 


+  5(3  0cc , 


(5) 


LAT  dt 

where  n(t)=nco{[+<xt),  z  =  L(l+a  t)  and  a  <*10  5  adaa1  typically.  Besides,  we  suggest  that 

field  frequency  is  near  the  cut-off  frequency  of  the  fiber.  Then,  the  transverse  waveguide 
parameter  U  is  [3]: 


U  -  V  exp(--!-)  «  2,4exp(--^-) . 
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At  the  same  time,  (3  =  knc 
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and  the  propagation  constant  of  the  FFEX ,  mode  is 


(30  »  knco  (l  -a  t) 


1  -  2A(2,4): 
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where  V  ~  V (l  +a  t ) . 

On  the  one  hand,  the  CV  vortex  has  U ^  ~  V  and 

Pcv  »P  =knJl-2A), 


(8) 


or 


5(30=Ho(2A)  l-(2,4): 
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Consequently,  the  temperature  sensitivity  is  defined  as 
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Acp  ,  ,  aL  J V-l\v-2) 
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On  the  other  hand,  it  can  be  considered  that  2a  kA  o=  1  and  _^P_  ~  4 — j-or  y_3 

L AT  i  adaa 

It  means  that  a  heating  of  one-meter  fiber  by  1°  C  rotates  the  intensity  zero  on  1,  27  turn  around 
an  optical  axis. 

2)  All  three  modes  propagate  along  a  parabolic-index  fiber.  Then  A(3  =  0  [1,3]  and  thetrajectory 
ofzero  is 

F  F 

x- - cosfLz,  y= - sinfLz.  (11) 

G+ 1  G-l 

F  F 

Thus,  the  vortex  trajectory  has  an  elliptic  form  with  semiaxes  a  = - ,  b  = - .  The 

G  +  l  G-l 

temperature  sensitivity  does  not  change  compared  with  the  first  case. 

3)  All  three  modes  propagate  along  a  step-index  fiber.  The  polarization  correction  of  the  CV  is 
[3] 

_  J2A)*  WU^Kfiv) 

Mcv  p  ]/3  K0(W) 
where  Kn  (.v)  is  the  McDonald  function. 

In  order  to  take  into  account  the  temperature  changes  of  the  polarization  correction,  one 
assumes  that  U  ~  U id  +  {IA)^  AV ,  where  AV  the  detuning  of  U  relative  to  the  cut-off 

frequency.  Besides,  K0 {W ) «  -In  VK  +  (in  2  - y ),  y  «  0,58  ,  K,  (w) «  1  /W  and  W  «  -J2VAV  . 
Whence, 

(2A)^2  1 

613^-2^ - — - - -.  (13) 

pV  ln(2VAy) 

and 

AO_  (2A  )3AWU2K1(W) 


AP=-- 


2p  y3  k2(w) 


But  K2(w)~  — — ,  consequently  A(3  «lf2  — >  0 . 


Assuming  that  Fl  =  Re xp(~]  / 2VR2>  e  F  =  exp(-l/2 VR2)  one  finds  the  trajectory’s 


equation: 


x  =  -7 r-^r— rcos5^z’  Y  =  7T— 7^— rsin8(30z 
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(a) 


(b) 


Fig.l.  The  trajectories  of  OV  in  the  temperature  sensor  at:  (a)  F=l,  G=3,  At  =  1  °C  , 
z  =100  m.  The  dependence  of  the  angular  position  cp  on  temperature  at  F=l,  G=0. 1 , 
z  =  0.3  m  (b). 


The  variable  amplitudes  in  eq.(15)  point  to  a  precession  movement  of  the  OV  trajectory. 
Provided  8(3  o  »  8(3<r  and  G«l,  the  temperature  sensitivity  will  be  really  the  same  as  in  the 
first  and  second  cases.  The  Fig.l  and  Fig.2  show  the  OVs  trajectories  and  the  intensity 
distribution  as  function  of  temperature. 


The  s  temperature  sensitivity  defined  from  Fig.  lb  is  about 


A(p  rad 
LAT ~  '  Km 


and  has  a 


good  agreement  with  the  obtained  result  for  the  first  case. 

The  experiment  is  brought  about  by  means  of  the  set-up  with  the  single  interferometric 
arm.  The  input  and  output  ends  of  the  fiber  are  positioned  on  the  three-coordinate  tables.  The 
laser  spot  at  the  plane  of  the  fiber  input  was  controlled  up  within  0.5  pm  and  observed  by  the 
reflected  light  microscope. 


(a) 


30 


(C) 


(d) 


Fig.2.  The  intensity  distribution  of  the  mode  field  viewed  at  the  different  temperature  values: 
25°  C  (a,c)  and  26°  C  (b,d),  respectively:  G  =0.2,  F  =1  (a,b);  G  =1,  F  =1  (c,d) 


The  low-birefringent  fiber  was  made  by  CVD  method  and  consists  from  G  e02 1  Si02 
core,  the  B20  /  Si 02  cladding  and  the  silica  substrate.  It  had  the  circular  geometry  with  the  core 

radius  P  =3.4 pm.  The  linear  birefringence  was  less  than  10“6.  This  fiber  was  designed  to  be 
single-mode  at  X  =  0.85  pm  but  the  experiments  were  performed  at  0 .63  p m  wavelength  region 
where  the  fiber  guides  three  waves:  HEn  mode,  CVQ1  and  IV01  vortices.  The  length  30-m 
section  of  the  fiber  was  looped  in  several  loose  coils  and  was  placed  in  a  temperature  controlled 

water  bath.  The  output  fiber  end  was  magnified  by  the  20x  microobjective,  recorded  by  the 
CCD-matrix  with  512x582  sensible  elements  and  viewed  on  the  TV  monitor  to  be  processed  by 
the  computer  -  Pentium  3  (800).  The  angular  positions  of  the  zero  intensity  in  the  field  image  was 
measured  as  a  function  of  temperature 


Fig.3.  The  experimental  dependence  of  the  rotation  angle  (p  of  the  OV  center  on  the  temperature 


t 0  C  .  The  experimental  accuracy  is  about  1.5%,  the  fiber  length  is  30  m. 


The  circular  analyzer  placed  after  the  fiber  output  suppresss  partially  the  «parasitic»  noise 
inserted  by  the  IV  mode.  The  value  of  the  noise  suppression  was  about  20  dB. 

The  Fig.3  demonstrates  the  experimental  curve  of  the  temperature  sensitivity.  The 
sensitivity  of  the  device  can  be  seen  to  be  a  linear  function  of  the  temperature  at  least  at  the  range 

rad 

from  20  to  40  °C  the  measured  temperature  sensitivity  is  about  1.03 -  and  in  good  agreement 

Km 

with  the  represented  above  theoretical  estimations.  The  reproducibility  of  the  experimental  results 
for  different  pieces  of  the  given  fiber  is  about  3.4%. 
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4.3  Optical-vortex  pulses  in  OVFS 
4.3.1  General  analysis 


Behavior  of  the  OV  pulses  in  a  fiber  inserted  to  an  interferometer  is  the  subject  of  the 
peculiar  attention  due  to  special  properties  of  the  interferometric  sensor.  At  the  same  time  this 
problem  has  not  been  yet  studied. 

The  aim  of  this  paragraph  is  to  give  the  theoretical  estimation  and  computer  simulation  of 
a  pulse  OV  propagation  in  an  optical  fiber  without  perturbations. 

Consider  a  quasimono  chromatic  pulse  in  the  form 


E(t,r)=E0exp \--j  exp(-  /  co0  t)F(r). 


T  2 

V  'o  J 


(1) 


Here  F  (r)  is  a  profile  function  of  the  beam,  r  =  *Jx2  +  y2  ,  T0  is  the  pulse  duration,  E  0  is  the 
amplitude,  CO0  is  the  central  frequency. 

The  spectral  decomposition  of  the  beam  gives 


E  (co )  =  |  exp  (i  co  t  )E  [t  ')dt ' , 
or,  taking  into  account  (1) 


(2) 
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where  V0  =  —  «  CO0 . 

' o 

Thus,  the  spectral  distribution  has  a  form  of  a  narrow  peak  in  the  vicinity  of  CO 0  with  the 
bandwidth  v0. 

Let  us  take  that  the  pulse  comes  to  the  fiber  input.  Naturally,  it  excites  the  whole  mode 
spectrum.  If  the  only  guided  OV  is  excited  in  the  fiber  then  the  wave  field  propagating  along  the 
fiber  is 


E  -C(m)E(t,  r,z)  exp(?P0)z-/cof), 
where  p0)  =  P0)(oj)  and 
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.  The  fiber  field  at  the  point  z  at  moment  t  equals 

©<p (/ (3  (co ) z — / co t) F  (r)  . 


(5) 


Because  the  interval  is  very  narrow,  the  range  2  V0  is  cut  out  from  the  all  integration  area.  Since, 
we  can  write 


32 


(3  (to)  «  (3  (cd0 )  + 13 ' (co0 ) (cd  -co0)+-^P"(co  -a>0)2 
After  the  integration  of  (5)  by  the  Laplace  method  we  obtain 


E(t,r,z ) 
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dp  (co) I 


.  Here  \/  is  the  group  velocity  of 


the  pulse  and  the  constant  D  characterizes  the  pulse  broadening  [1]  in  the  fiber. 

In  order  to  estimate  the  pulse  broadening  and  the  group  velocity  of  a  fiber  OV  let 

8 1 .  =2  —  DHASoj  .  However  8  t.  =  — — y  8(0  since  we  find  D  =  —  ^  ^ —  DV  A, 

J  cc o  1  lvcfo)2J  2  dco2  cco 

where  nco  is  the  refractive  index  of  the  fiber  core,  A  is  the  height  of  the  refractive  index  profile, 

V  is  the  waveguide  parameter:  V2  =  U  2  + 1/1/ 2 ,  U  and  W  are  the  core-  and  the  clad-  waveguide 
parameters,  respectively.  Besides,  the  group  velocity  is  evaluated  for  a  step-index  fiber  as  [2] 


C  U2 

V — i  1  + Arrr  1-2 


and  /  >  1,  K/(W)  is  the  McDonald’s  function. 

For  the  typical  low-order  CV  in  the  step-index  fiber  with  V=5,  we  have  D=-0.01 ,  U=3. 1 5,  W=3.9, 

r2 

T0  =10“8c,  n  =148,  n  ,  =1.46,  A=  0.014,  A,  =  0.63  pm,  whence  D  «116xl0“27 — j, 

m 

TYl  Z 

v  ~  2.07  xlO8 —  and  LD  ~  8.6x10 10m .  Consequently,  not  very  short  pulses  have  —  «  0  . 


At  the  OVFS  output  the  wave  field  is  the  sum  of  the  signal  Esand  the  reference  ER 


beams: 


Es  ocE0exp  - 


t-  — 


F  fr)exp(-/co  t)’E«  -g0exp-V  ^  ~  E  (r)exp(-/cot),  (9) 


where  Zs ,  ZR  are  the  lengths  of  the  signal  and  the  reference  arms. 
The  total  field  intensity  at  an  observation  plane  is  written  as 
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where  cp  is  the  azimuth  angle,  A((3  z)  =  (30ZS  -  k  ZR  . 


As  it  can  be  seen  from  (10),  two  processes  affect  the  pulses  and  the  interference  pattern. 
The  first  process  is  defined  by  difference  of  time  of  the  pulses  propagation  through  the 
interferometer  arms.  The  second  process  is  the  frequency  dispersion.  As  the  result,  the  image 
contrast  of  the  interference  patterns  diminishes. 

To  evaluate  the  image  contrast  let  us  pass  to  dimentionless  variables: 

*  =  7 WT.  zo  =  vgTo  =  2.7m,  ?„=V2y-,  = 

W  Z  /  o  Zq  Zq 

It  should  be  taken  into  account  that  the  photodetector  records  time-avereged  value.  In  our  case, 
the  pulse  duration  T0  is  time  of  averaging.  Consequently,  the  observed  intensity  distribution  will 
be: 

</>=(/s)  +  </«)  +  (/«s>C0^A(pz)-/q>],  (11) 

where  (/s )  =  J^{erf  (V2  - ^ )  +  erf  fe  )}f (r) , 

(/««)  =  if{erf  +erf  |fe+Cs)|exp-|fc-Cs)2  F(r). 

Thus  the  image  contrast  is 


K  = 


max  f  min 
max  J  min 


M 

(».)+(',) 


(12) 


where  /max  and  /min  are  the  maximum  and  minimum  intensities  of  the  interference  image 
respectively. 

The  Fig.l  illustrates  the  evolution  of  the  intensity  distribution  of  caused  by  interference 
oftwo  OV-pulses  transmited  through  the  two-arm  interferometer.  (The  patterns  of  the  Fig.l  are 
depicted  in  accordance  with  equations  (10),  (11)).  The  interferometric  sensor  is  tuned  to  the  zero- 
order  band  when  the  only  dark  straight  line  is  observed.  The  angle  direction  of  the  line  is 
proportional  to  the  phase  shift  between  the  beams. 
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(a)  (b) 


Fig.2.  The  dependence  of  the  image  contrast  K  on  the  dimensionless  length  and  ^R  of  the 
interferometer  arms. 


As  the  difference  between  the  beam’s  optical  paths  increases,  the  image  clarity  sharply 
reduces  and  practically  vanishes  at  =  2.5,  C,R  =  1.  In  fact,  the  interference  band  is  spread.  It 


corresponds  to  the  following  optical  path  ratio 


2.5.  The  clarity  increases  with  increasing 


the  pulse  duration  T0 . 

It  is  convenient  to  characterized  the  clarity  by  means  of  the  image  contrast  defined  by  Eq. 
(12).The  Fig.2  represents  the  set  of  curves  describing  the  dependency  of  the  image  contrast  on  the 
dimensionless  and  ^R  parameters.  Negative  or  positive  -values  corresponds  to  advancing 
or  delaying  of  pulses  relative  to  one  another.  The  behavior  of  the  contrast  curves  shows  that 
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increasing  the  phase  shift  between  the  beams  and  decreasing  the  pulse  duration  causes  the 
spreading  of  the  interference  image.  The  sharp  peaks  on  the  curves  at  the  small  phase  differences 
are  gradually  reduced  at  the  expense  of  either  the  optical  path  increasing  or  the  pulse  duration 
decreasing  and  they  can  not  essentially  effect  the  interference  pattern  quality. 

Thus,  employment  of  short-laser  vortex-pulses  in  two-arm  interferometric  sensors  is  not 
efficient.  Indeed,  there  exist  two  main  optical  processes,  which  destroy  an  interference  pattern  in 
the  OVFS.  The  first  process  is  affected  by  advancement  or  retardation  of  pulses  relative  to  one 
another  so  that  their  contours  are  partially  overlapped  or  are  not  overlapped  at  all.  The  second 
process  is  associated  with  the  broad  spectrum  of  monochromatic  waves  that  corresponds  to  any 
short  pulse.  Both  these  reasons  cause  the  scrambling  of  isolated  interference  patterns,  and 
aggravation  of  the  image  quality. 
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4.3.2  Super  sensitive  one-arm  OVFS 


In  the  work  [1]  the  problem  on  a  super  sensitive  temperature  measurement  by  means  of 
an  one-arm  fiber-optic  interferometer  on  the  base  of  the  fundamental  HEn  eigen  mode  of  a 

single  mode  fiber  has  been  discussed.  That  temperature  sensor  represents  a  system  consisting  of  a 
single  mode  fiber  at  the  ends  of  which  two  special  mode  couplers  are  fitted.  The  only  linear 
polarized  light  component  passes  through  the  single  mode  fiber  if  an  electrical  vector  of  an  input 
beam  is  oriented  along  a  main  fiber  axis.  The  first  coupler  divides  the  initial  linearly  polarized 
mode  between  the  two  linearly  polarized  ones.  The  second  coupler  combines  these  modes 
yielding  the  mode  polarization  interference.  The  fiber  excitation  is  brought  along  by  a  tungsten 
lamp  with  a  broadband  of  a  light  spectrum.  Light  emitted  from  the  fiber  output  is  analyzed  by  a 
spectrograph.  As  the  authors  point  out  this  fiber  sensor  has  a  temperature  resolution  of 

3x  10“3  °C  for  a  wide  temperature  range. 

The  first,  it  should  be  noticed  that  the  above  mentioned  sensor  is  unstable  one  to  any 
slight  perturbation.  Indeed  the  single  mode  fiber  in  this  device  is  excited  by  a  broadband  light 
source  (a  tungsten  lamp).  But  a  single  mode  fiber  is  known  to  have  a  cut-off  frequency  for  high- 

order  modes  at  V  =  2.4,  where  1/  =  -j—  p-Jri.,  —  H;;  -  the  waveguide  parameter,  A,  is  a 


wavelength,  p  is  a  fiber  core  radius,  nco,nd  -  refractive  indices  of  a  core  and  a  clad, 
respectively.  Besides,  a  light  beam  with  a  wavelength  less  then  the  cut-off  wavelength  excites 
high-order  modes  in  the  fiber.  The  authors  make  use  of  the  F-SPV  birefringence  optical  fibers  that 
have  the  critical  wavelength  X  ~  0.85|lm .  In  these  fibers  the  H  En  H  E21  even  and  odd,  TE01 

and  TM  01  modes  are  realized  at  the  wavelength  X  —  0.74pm  at  which  the  authors  worked.  Any 
slight  perturbation  of  the  fiber  means  as  causing  an  energy  mode  conversion  that  changes  the  light 
spectrum  and  to  induce  a  mode  noise. 

The  second,  the  mode  polarization  couplers  at  the  fiber  ends  are  deformed  due  to  boiling 
or  cooling  and  change  the  mode  transformation  value.  It  causes  a  mode  noise  as  well. 

The  aim  of  the  given  section  is  to  study  the  new  fiber-optic  sensor  on  the  base  of  an 
optical  vortex  devoid  of  the  above  mentioned  deficiencies. 

The  main  idea  is  as  follows.  If  the  optical  fiber  excites  the  light  beam  with  a  topological 
charge  /  and  a  circular  polarization  a  having  different  signs,  there  are  bom  three  mode  types  in 


the  fiber,  namely  a  steady  CV  vortex,  an  unstable  IV  vortex  and  a  fundamental  HEf°}  mode.  The 


CV  vortex  and  HE{°} 


mode  have  the  same  polarization  state  as  the  initial  light  beam  while  the 
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polarization  state  of  the  IV  vortex  may  be  changed.  The  circular  polarization  analyzer  separated  a 
circular  polarized  light  gets  to  fit  at  the  fiber  output.  This  analyzer  suppresses  one  of  circular 
polarizations.  Let  the  analyzer  will  be  tuned  in  the  IV  mode  polarization.  Then  the  CV  and  HE 
modes  are  suppressed  and  the  IV  vortex  turns  out  to  be  the  work  wave.  As  the  matter  of  fact  that 
the  IV  vortex  is  a  sum  of  symmetric  TE01  and  TM  01  modes  transmitting  with  different  phase  and 
group  velocities.  In  fact,  we  are  concerned  with  two  arms  interferometer  but  its  two  arms  are 
combined  into  one  directed  line  -  an  optical  fiber.  It  should  be  worth  noted  that  our  method  in 
contrast  to  the  previous  one  is  not  needed  the  special  birefringence  fiber  types  but  it  can  be  used 
an  ordinary  commercial  optical  fiber. 

Consider  this  method  in  detail.  The  wave  function  of  the  light  field  after  passing  through 
the  fiber  and  the  circular  polarization  analyzer  is 


VFH  =  ¥0U)exp(-q>)sinf^^-z  ]  exp(/p~z)F  (r), 


(1) 


where  the  optical  fiber  is  suggested  to  be  excited  by  a  right-hand  circularly  polarized  optical 
vortex  with  the  topological  charge  1=1 ;  the  analyzer  cuts  off  the  right-hand  circular  polarization 
while  the  left-hand  polarized  light  component  passes  through  it  without  losses.  The  sign  (-)  means 


a  left-hand  circular  polarization,  A  (3  = 


^Ptm  ^Pte 


P  is  a  scalar  propagation  constant, 


8pTM  and  8(3^  are  polarization  corrections  to  the  scalar  propagation  constant  (3  ,  F  (r)  is  a  field 
profile  function,  *P0(  A,)  is  the  input  intensity;  the  light  field  propagates  along  z-axis. 


Fig.l.  The  spectral  response  S(X)  of  the  optical  fiber  with  IV  vortex,  A fT0  ~  10  9 ,  z  =1  m. 


In  a  general  case  the  difference  of  the  polarization  correction  for  the  TM  and  TE  modes 
may  be  written  in  terms  of  the  effective  refractive  indices  Aft  =  flTM  -  flTf  in  the  form 


AP  =  -^A/?U) 

A 


(2). 


The  dependency  A/T  on  a  wavelength  k  is  suggested  as  being  far  off  absorption  bands  and  to  be 
described  by  the  phenomenon  of  a  normal  material  dispersion. 

The  spectral  response  of  the  optical  fiber  is 


SU)  = 


Iwl3 


=  cos 


r-An(X)z  F2(r). 


(3) 


Assuming  the  normal  dispersion  has  a  weak  response  in  the  given  spectral  range  Aff[X)  ~  AfT0  it 
is  found  out  the  plotting  dependency  shown  by  Fig.  1 

The  curve  on  the  Fig.l  has  an  oscillation  form  and  equals  zero  at  the  points 

^-AfT(X), 

for  the  wavelength  A,  and  A,  respectively,  M  is  the  number  of  cycles  between  \  and  X2 .  N  is 
integer.  Excluding  z  from  (4)  we  have 
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(5) 
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For  the  neighboring  wavelengths  we  have  M=  1  and 


N 


K 


l 


(6) 


4 

The  temperature  control  is  brought  along  by  means  of  the  spectral  measurement,  i.e.  the 
temperature  changes  causes  the  wavelength  shift  and  are  recorded  by  a  spectrograph. 

dX 

Evaluate  the  spectral  temperature  sensitivity  — — .  From  (4a)  we  have 
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A  boron-doped  fiber  has 


dA  n 

aT 
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,  where  Ts  is  a  softening  temperature  so  as  for 


-r  -r  dz  _ _  ,  ._7  An  z 

room  temperature  T  —  Ts  =  800  C  and  ~z5.5x  10  [2],  If  M=1  we  have  - ^"~^i 


2  N  +■ 


and  the  spectral  temperature  sensitivity  is 

dX 


dT 


yT-Ts 


+  55x10 


i-7 


(8) 


or 


at  wavelength  X:  —  0.63| im. . 


On  the  base  of  the  above  carried  over  consideration  the  follow  conclusions  may  be  drown. 


The  spectral  temperature  sensitivity  of  the  one-arms  fiber-optic  interferometer  does  not 
depend  explicitly  on  a  fiber  length.  It  is  necessary  only  the  temperature  changes  must  cause 
changes  of  an  optical  fiber  length  no  less  than  one  wavelength. 

The  advantages  of  the  given  method  are  follows. 

At  first,  when  being  a  sum  of  the  symmetric  TM  and  TE  modes  the  IV  vortex  is  not  very 
sensible  to  changes  of  a  cut-off  frequency  at  least  as  compared  with  H  En  mode.  Indeed,  as  far  as 


this  fiber-optic  system  maintains  only  IV  vortices  it  is  responded  to  birth  and  death  events  of  the 
TM  and  TE  modes  exclusively.  The  first  mode  pair  is  bom  at  V=2.405  while  the  second  one  does 
it  at  V=5.52.  This  waveguide  parameter  corresponds  to  the  cut-off  wavelength  Xc  =0.347  mem . 
It  means  that  this  wavelength  is  not  analyzed  by  the  experimental  set-up.  This  thing  extends  the 
dynamics  range  of  the  device  and  reduces  the  measurement  error. 

The  second,  our  experimental  set-up  has  not  the  spiral  polarization  couplers  since  the 
measurement  errors  due  to  a  temperature  response  of  the  couples  are  vanished. 


The  high  spectral  sensitivity  about 


associated  with  the  temperature 


resolution  up  to  2..4x  10“3  °c  enables  us  to  originate  super  sensitive  arrangement  of 
temperature  and  other  thermodynamics  values. 
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4.4  New  perspective  approaches  and  their  estimations 

It  was  shown  above  that  low-mode  OVFS  possess  the  greatest  resolution  among  optical 
fiber  sensors,.  The  basic  mechanism  is  the  temperature-induced  changes  of  the  phase  and 
polarization  states  of  a  OV  beam  guided  through  thenfiber.  Unfortunately,  the  transmission  of  OV 
through  the  conventional  commercial  fiber  is  very  sensitive  to  external  mechanical  perturbations, 
so  that  the  OV  is  decayed  into  different  eigen-modes  and  a  mode  noise  is  induced  in  the  device. 

It  has  been  well  known  technique  of  fabrication  of  the  single-mode  fiber  at  the  expense  of 
strong  stresses  inserting  into  a  fiber  material.  Such  optical  fibers  can  preserve  the  mode 
polarization  even  if  greatest  ambient  forces  affect  them. 

In  the  given  paragraph  we  shall  consider  some  fiber  constructions  enable  to  preserve  a 
single  OV  at  large  fiber  lengths,  which  are  designed  on  the  base  of  single  polarization-preserving 
fibers. 

Assume  that  there  are  two  kinds  of  optical  anisotropy  in  some  optical  material,  namely  a 
linear  birefringence  A nL  and  a  circular  one  Aflc  ,  where  An  is  a  difference  between  refractive 
indices  of  eigen-  modes.  In  terms  of  plane  waves,  the  response  of  that  complex  birefringent 
medium  is  expressed  through  the  effective  difference  of  the  refractive  indices  A nR  in  a  form 
[1,2]: 

A  nR  =  ^Anf+Anf  (l) 

Respectively,  if  the  fiber  has  a  complex  birefringence,  the  polarization  corrections  to  propagation 
constants  of  eigen-modes  appear.  If  A nL  » A flc ,  the  principal  process  is  the  linear 
birefringence  and  vice  versa  at  AflL  «  A flc  . 

The  main  noise  in  the  OVFS  gives  the  IV  mode.  Its  polarization  correction  is  associated 
with  the  linear  topological  birefringence  8|3L  [3],  while  the  correction  of  A(3C  for  the  CV  mode 
is  unaffected  by  noise  phenomena,  whichoccur  due  to  the  circular  birefringence.  Consequently,  if 

the  strong  circular  birefringence  takes  place  in  the  fiber  (8|3C  »  8(3L ),  the  mode  noise  in  the 
fiber  would  vanish. 

It  has  been  also  known  [4]  that  the  so-called  “spun  fibers”  preserve  the  circular 
polarization  of  the  fundamental  H  En  mode.  They  have  an  elliptically  screwed  core  and  are 

fabricated  by  spinning  of  the  platform  during  fiber  drawing  to  impact  the  permanent  twist  of  a 
fiber  and  thus  to  restore  in  average  the  circular  symmetry  of  the  waveguide  structure.  However, 
the  transformation  processes  of  guided  OVs  in  such  fibers  are  not  yet  studied. 

In  order  to  estimate  the  main  features  of  the  OV  transformation  process  we  have 
considered  a  fiber  with  a  spirally  screwed-core  (see  Fig.l). 


(b) 


39 


Fig.l.  The  spirally  screwed-core  fiber:  (a)  the  single  spiral;  (b)  the  dual  spiral. 


A  beam  transmited  through  the  spiral  trajectory  gains  the  topological  Berry’s  phase  [5]. 
In  our  case  the  Berry’s  phase  per  one  turn  of  the  spiral  has  a  form: 


y  =  27t(l-cosf>)  =  271 
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where  r  isthe  radius  of  the  spiral,  h  is  the  spiral  pitch.  The  specific  topological  phase  per  unit 
length  is 
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The  magnitude  of  the  specific  topological  phase  characterizes  the  polarization  correction  to  the 
scalar  propagation  constant  [6] :  I  =  8(3C  .  Evidently,  for  the  case  8(1-  »  8(1  ,  the  linear 
topological  birefringence  will  be  suppressed  in  the  fiber,  and  the  mode  noise  vanishes.  For 

simplicity,  we  assume  that  |S(3j  ~  |S(3C|  in  a  step-index  fiber  [7]  and 


8|3C 


V2A 
pV  ' 
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where  p  is  the  fiber  radius,  A  is  the  height  of  the  refractive  index  profile,  V  is  the  waveguide 
parameter  and  azimuth  index  1=1. 


Fig.2.  The  characteristic  curve  of  the  spirally  screwed-core  fiber.  The  scale  r  and  h  is  in 


micrometers.(a).  The  same  characteristic  curve  in  the  coordinates  u  =  10 
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The  Fig.2a,b  represents  the  graphs  of  the  characteristic  curves. 
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The  first  curve  in  the  Fig.2a  was  obtained  from  the  inequality  8[3C  »  8(3;  and  eq.  (3),  (4)  for 

typical  fiber  parameters:  A  =  10  2 ,  p  =  4|lm  ,  V=3.5.  It  is  associated  with  the  parameters  of  the 

spiral  core  only,  while  the  second  curve  of  the  Fig.2b  depends  both  from  the  parameters  of  the 
straight  fiber  itself  and  from  its  core  torsion.  The  values  of  the  spiral  radius  r  and  its  pitch  h  above 
the  curve  of  the  Fig.2a  correspond  to  the  parameters  of  the  spiral-core  fiber  preserving  the  guided 
OVs.  At  the  same  time,  the  values  settled  above  the  curve  of  the  Fig.2b  point  to  the  admissible 
values  of  the  straight  fiber  parameters  at  the  given  values  of  the  torsion. 

As  a  matter  of  fact,  the  real  physical  processes  in  the  fiber  are  more  complicated  than 
presented  above.  Indeed,  the  spiral  fiber  is  the  waveguide  with  an  imperfection  of  the  refractive 
index  along  the  optical  axis.  There  is  longitudinal  modulation  of  the  refractive  index  in  the  fiber, 
and  hence  there  arises  a  counter  waves  due  to  the  Bragg  scattering.  Moreover,  azimuthally 
symmetric  TE  and  TM  modes  are  not  realized  in  this  system  because  of  violation  of  waveguide 

symmetry.  Instead  of  them  it  appears  the  new  C  V  s  with  a  slightly  different  polarization 
correction  than  the  CV  has.  As  a  result,  any  slight  thermodynamical  perturbations  of  temperature, 
pressure  and  so  on,  the  spiral  pitch  changes  cause  violation  of  the  Bragg’s  condition  and  detuning 
of  the  polarization  corrections  of  the  CV  and  C  V  s  vortices.  If  a  polarization  analyzer  is 
positioned  after  the  fiber  output,  there  arise  strong  variations  both  of  light  intensity  and 
thespectrum  at  the  output  of  the  OVFS. 

We  found  also  that  the  similar  properties  are  to  the  elliptically  screwed -core 
fiber  (see  Fig.3) 


Fig.3.  The  elliptically-screwed  core  fiber. 

It  must  be  emphasized  that  the  described  above  OV  traits  are  mainly  associated  with 
phase  properties  of  a  wave  than  features  of  the  OV  beams.  The  point  is  that  that  guided  OV  in  the 
fiber  is  itself  an  optical  gyroscope  similar  to  the  mechanical  top  in  a  frame  in  contrast  to  the  OV 
in  free  space  that  is  similar  to  the  top  without  a  frame  [8].  The  main  distinctive  trait  of  the  guided 

OV  is  the  polarization  correction  8|3  to  the  scalar  part  of  the  propagation  constant  (3  .  The 
polarization  correction  8(3  has  been  shown  [8]  to  be  immediately  associated  with  the  spin-orbit 
interaction  as  the  same  as  that  takes  place  for  the  mechanical  gyroscope  and  other  mechanical 
devices  of  the  inertial  navigation.  Any  variations  of  an  angular  position  of  the  fiber  coil  stimulate 
changes  of  the  polarization  correction.  Moreover,  those  changes  depend  both  on  a  topological 

charge  and  a  rotation  rate.  The  magnitude  of  the  polarization  correction  8|3(G)is  proportional  to 
the  area  enclosed  by  a  single  coil  of  a  fiber  during  rotation  and  to  the  coil  number.  If  two  OVs 
formed  as  a  topological  dipole  propagate  along  the  rotating  fiber,  their  centers  begin  to  roll 
around  the  optical  axis  also.  The  rotation  rate  of  the  fiber  coil  is  proportional  to  the  rotation  rate 
of  the  topological  dipole.  In  contrast  to  the  conventional  fiber-optical  gyroscopes  in  configuration 
of  the  Sagnac  interferometer  proposed  OV  device  eliminates  the  counter-propagated  waves  and, 
consequently,  eliminates  all  experimental  errors  inherent  to  a  such  OV  interferometer. 
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4.5  Conclusions  to  chapter  4 


We  have  made  the  analytical  review  and  have  shown  the  advantages  and  imperfections  of 
optical  sensors  on  the  base  of  both  single-mode  and  low-mode  fibers.  As  a  rule,  a  high  sensitivity 
of  these  devices  is  immediately  combined  with  errors  at  the  expense  of  intensity  fluctuations  of  a 
laser  source.  Besides,  experimental  results  obtained  by  these  sensors  depend  linearly  on  the  beam 
phase  difference  only  at  narrow  range  of  magnitudes. 

We  have  developed  and  tested  the  specific  physical  approach  to  creation  of  the  novel  type 
of  optical  sensor:  the  OVFS,  lacking  a  number  of  imperfections  of  the  mentioned  above  schemes. 
This  approach  considers  the  OV  interferometer  as  the  complete  wave  system  involving  wave 
transmitting  processes  through  the  OV-generated  unit,  the  vortex-guided  unit,  the  vortex- 
recovered  unit  and  on  the  final  stage  an  interference  image  is  subjected  to  a  computer  analysis  by 
the  processing  unit.  In  accordance  with  our  approach,  the  OV  is  employed  as  the  signal  or 
reference  beam.  The  image  obtained  as  a  result  of  the  waves  superposition  has  the  form  of  the 
interference  spiral  rotating  around  its  axis  due  to  the  phase  difference  between  the  OV  and  a 
Gaussian  beam.  The  computer  processing  of  the  spiral  image  instead  of  the  measurement  of  the 
total  intensity  enables  to  avoid  influence  of  intensity  fluctuations.  An  angular  position  of  the 
spiral  image  is  proportional  to  the  phase  difference  at  a  broad  range  of  magnitudes.  Two  possible 
interferometric  schemes  were  considered:  (i)  with  an  OV  as  the  signal  wave  guided  by  the  low¬ 
mode  fiber  and  (ii)  with  an  OV  as  the  reference  wave  transmitted  through  free  space  and  the 
single-mode  fiber  as  the  signal  arm.  The  former  device  has  temperature  sensitivity  about  180 
rad  rad 

-  while  the  latter  one  about  110  - .  The  specific  technique  of  the  image  processing  has 

Km  Km 

been  operated  so  that  the  measurement  error  of  the  device  was  only  0.01  rad  and  less  and  the 
reproducibility  was  up  to  2  %.  However,  the  former  interferometer  respondes  rather  strongly  to 
ambient  mechanical  perturbations  and  the  influence  of  coupled  OVs  implemented  in  the  fiber. 
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This  mode  mixture  was  experimentally  and  theoretically  revealed.  The  coupled  OVs  are 
suppressed  and  the  single  guided  OV  is  recovered  almost  without  any  energy  losses  by  the  two- 
lens  system.  Also  we  have  developed  the  technique  of  OV  generation  by  an  optical  wedge. 

Transmission  through  a  fiber  of  short  laser  pulses  is  restricted  by  pulse  duration.  On  the 
base  of  both  the  analytical  consideration  and  computer  simulation,  we  found  that  the  contrast  of 
an  interference  pattern  is  characterized  not  only  by  pulses  mutual  correlation  in  space  and  time 
but  also  by  pulse  dispersion.  Two-arm  OVFS  cannot  operate  as  optical  sensors  in  the  pulse 
regime.  Nevertheless,  we  have  developed  and  tested  the  one-arm  OVFS  suitable  both  for  cw  and 
pulse  regimes.  In  this  scheme,  the  CV,  and  H  En  modes  propagate  along  the  only  fiber.  The  IV  is 

partially  suppressed  by  the  circular  polarization  analyzer.  The  CV  and  HEn  modes  have 
different  responses  to  temperature  influence  and  thus  represent  two  interferometric  “arms” 
launched  into  the  single  fiber.  The  angular  position  of  the  zero  point  of  intensity  characterizes  the 

rad 

value  of  temperature.  Such  one-arm  OVFS  has  low  temperature  sensibility,  nothing  but  2 - . 

Km 

However,  we  have  theoretically  developed  a  specific  spectral  technique  raising  the  temperature 

•  •  •  rad 

sensitivity  up  to  200 - and  more. 

Km 

We  have  studied  the  physical  mechanism  of  OV  decay  in  the  conventional  fibers  and 
offered  two  novel  fiber  types  preserving  OVs.,  namely  fiber  with  the  spirally  screwed  core  and 
the  elliptical  fiber  with  the  twisting  core,  or  «a  spun  fiber».  Our  estimations  show  that  OV  sensor 
on  the  base  of  the  OV-preserving  fiber  is  not  merely  high  sensitive  interferometer  of 
thermodynamical  values  but  it  can  be  considered  as  the  sensor  for  inertial  navigation. 

Conclusions  to  the  project  P  051 


During  one-year  Partner  Project  P-051  a  set  of  essential  contributions  were  made  into  the 
understanding  of  the  nature  of  optical  vortices  in  free  space,  optical  fibers  and  nonlinear  media,  as 
well  as  optical  vortices  assessment.  Most  important  results  are  listed  below. 

Optical  vortices  understanding 

1.  New  interference  phenomena  with  optical  vortices. 

-The  structure  of  transversal  optical  vortex  around  edge  dislocation  was  analyzed.  It  was 
shown  that  interference  of  2D  Gaussian  beams  produces  pair  of  edge  dislocations.  Their  collapse 
and  “unfolding”  is  possible  in  saddle-vortex  topological  reaction.  Very  important  for 
understanding  of  edge  dislocations  was  the  consideration  of  phase  velocity  variation  in  the 
vicinity  of  zero-amplitude  line.  It  was  established  that  the  transversal  vortex  rotates  toward  the 
region  with  higher  phase  velocity. 

-It  was  predicted  theoretically  and  realized  experimentally  the  spontaneous  (natural)  birth  of 
OVs  by  overlapping  of  mutually  coherent  non-coaxial  co-propagating  Gaussian  beams  with 
different  parameters. 

2.  Diffraction  and  self-restoration  of  spatially  truncated  OV  beams 

-It  was  shown  by  computer  simulation  and  experimental  verification  the  effect  of  self¬ 
restoration  of  OV  after  up  to  80  %  cross-section  screening  of  an  OV  beam,  including  its  zero- 
amplitude  axis,  during  free  propagation  behind  the  screen.  This  new  nontrivial  effect  is  a  result  of 
redistribution  of  the  residual  orbital  angular  momentum  in  the  beam  cross-section  due  to  its 
rotation  and  diffraction  behind  the  screen. 
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3.  New  type  of  topological  charge  and  the  orbital  angular  momentum  transformation 

-The  observed  transformation  of  OV  beam  passing  through  a  cylindrical  lens  in  the  focal 
plane  and  vicinity  was  found  to  be  quite  unusual  and  unexpected.  It  was  shown  that  orbital 
angular  momentum  of  incident  OV  beam  is  conserved  contrary  to  the  topological  charge  of  the 
axial  OV,  which  changes  its  sign  during  passing  through  the  focal  plane.  This  effect  needs  more 
deep  theoretical  consideration. 

-Theoretically  calculated  effect  of  the  topological  charge  inversion  during  crossing  of  an  OV 
“dipole”  and  circular  edge  dislocation  was  revealed. 

Optical  vortices  assessment 

1.  Mode  separator  for  a  beam  with  an  off-axis  optical  vortex 

-The  technique  for  spatial  separation  of  modes  with  even  and  odd  values  of  modal  index  Q 
was  proposed  and  realized  on  the  basis  on  two-arm  interferometer  with  additional  Gouy  phase 
shift  in  one  arm  equal  -  Q  5.  This  simple  technique  allows  to  separate  Laguerre-Gaussian  modes 
and  “purifying”  in  high  extent  an  OV  beam  from  parasitic  coherent  background. 

2.  Multichannel  optical  manipulator  for  microparticles 

-The  physical  base  for  two-channel  optical  manipulator  for  microparticles  has  been 
developed.  The  manipulator  was  "constructed"  as  a  focused  singular  beam  possessing  two  OVs, 
with  only  one  easy  operating  parameter  for  management  of  the  OVs  positions  within  the  beam. 
Elaborated  computer  program  simulates  the  process  of  manipulation. 

-A  new  holographic  technique  was  elaborated  for  creation  of  OV  beams  with  fractional 
topological  charges.  It  was  predicted  and  shown  by  computer  simulation  that  SHG  of  OV  beam 
with  fractional  topological  charge  on  the  fundamental  frequency  possesses  quite  new  features. 
Namely,  the  topological  charge  of  SH  of  such  OV  beam  depends  in  crucial  way  on  the  distance 
between  the  nonlinear  element  and  the  computer-synthesized  hologram.  This  effect  can  be  used 
for  management  of  SH  of  OV  beam. 

3.  The  new  high-resolution  technique  for  testing  of  OV  beam  wave  front 

-  A  new  technique  for  testing  of  OV  helicoidal  wave  front  in  most  delicate  angular  range 
inside  the  angle  of  diffraction  divergence  was  elaborated  and  realized.  It  is  based  on  multi-beams 
interference  of  OV  beam  put  into  a  leaky  waveguide.  This  technique  deserves  subsequent 
development  and  application. 

4.  Fiber-optical  sensors  based  on  optical  vortices 

The  OV  temperature  sensors  on  the  base  of  different  types  of  optical  fibers  have  been 
theoretically  and  experimentally  studied.  The  OV  fiber-optical  sensor  interferometer  consists  of 
the  main  units:  the  excitation  system,  the  vortex  guided  system,  the  light  output  system  and  the 
detecting  -  computer  processing  system.  The  operation  of  the  each  unit  was  analyzed.  Two  types 
of  interferometric  configuration  were  considered:  one-  and  two-arm  interferometers.  The 
advantages  and  disadvantages  of  these  schemes  were  analyzed  and  the  ways  of  their  improvement 
were  designed. 

5.  The  new  technique  for  OV  beam  generation 
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-The  technique  of  efficient  generation  of  optical  vortices  in  free  space  was  developed. 
The  OV  beam  with  efficiency  up  to  70%  of  conversion  from  a  fundamental  Gaussian  beam  by 
diffraction  on  a  transparent  optical  wedge  was  obtained.  The  OV  generation  process  was 
considered  in  detail  both  experimentally  and  theoretically.  The  optical  and  mechanical 
arrangement  was  elaborated  for  the  optical  fiber  excitation. 

6.  Estimation  of  laser  short  pulses  propagation  in  OV  fiber  sensors 

-The  coherence  and  duration  of  laser  pulses  is  shown  to  limit  the  OV  interferometer  sensitivity. 
There  are  two  physical  phenomena  decreasing  the  sensitivity  of  the  vortex  devices.  The  first  one 
is  that  the  optical  path  difference  between  the  interferometer  arms  must  be  less  than  the  pulse 
length.  The  second  one  is  that  a  short  pulse  causes  the  spreading  of  an  interference  pattern  due  to 
a  broadened  spectral  range.  The  analysis  of  the  influence  of  these  limitations  was  performed  and 
computer  simulation  was  realized. 

7.  Experimental  realization  of  OV  fiber  sensors 

-The  one-arm  OV  fiber-optical  sensor  was  experimentally  developed.  The  temperature 

sensitivity  of  this  device  is  up  to  1  -  2  rac ^  .  The  two-arm  OV  fiber-optical  sensor  was 

Km 

experimentally  developed.  The  temperature  sensitivity  of  different  types  of  these  devices  is  up  to 

110-180  rad  and  the  experimental  error  is  up  to  3xlO~4K.  The  linearity  of  the  temperature 
K  -m 

measurement  and  reproducibility  was  verified.  The  test  protocols  of  the  examination  of  the 
devices  are  prepared. 
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